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Development  of  an  Airborne  Sea  Ice  Thickness  Measurement  System 

and  Field  Test  Results 

AUSTIN  KOVACS  AND  J.  SCOTT  HOLLADAY 


INTRODUCTION 

Airborne  electromagiretic  sounding  technology 
uses  multiple  sets  of  transmitter  and  receiver  coils 
that  function  as  dipole  antennas.  One  coil  of  the 
antenna  pair  is  moimted  forward  and  one  aft  on  a 
stiff  platform  securely  fixed  inside  a  rigid  cigar¬ 
shaped  shell  (bird).  Each  transmitter-receiver  coil 
pair  operates  at  a  single  frequency.  The  transmitted 
electromagnetic  field  induces  eddy  currents  in  a 
nearby  electrically  conductive  mass  such  as  sea¬ 
water.  As  a  result,  a  secondary  magnetic  field  aris¬ 
es  that  produces  an  analog  voltage  proportional  to 
the  secondary  magnetic  field  at  the  receiver.  The 
received  voltage  amplitude  and  phase  are  functions 
of  transmitter  coil  orientation,  position  with  respect 
to  the  receiver  coil  and  radiated  electromagnetic 
field  strength,  the  conductivity,  distance  and  geom¬ 
etry  of  the  conductive  surface  with  respect  to  the 
transmitter,  and  the  sensitivity,  orientation,  and 
distance  of  the  receiver  coil  with  respect  to  the  sea¬ 
water.  Since  the  electromagnetic  response  mea¬ 
sured  at  the  receiver  is  strongly  related  to  bird- 
seawater  distance,  an  accurate  measurement  of 
this  response  can  provide  a  very  good  estimate  of 
this  distance.  When  subtracted  from  the  bird  height 
above  the  sea  ice  surface,  this  distance,  as  measured 
with  a  laser  altimeter  mounted  in  the  bird,  gives  an 
estimate  of  ice  thickness. 

This  report  gives  an  accounting  of  the  effort 
made  to  implement  improvements  and  validate 
the  use  of  airborne  electromagnetic  induction  tech¬ 
nology  for  the  measurement  of  sea  ice  thickness. 
The  work  discussed  covers  the  period  from  June 
1986  to  February  1988  and  follows  the  pUot  study 
undertaken  in  1985  (Kovacs  et  al.  1987a).  In  the 
1985study,  airborne  electromagnetic  (AEM)  sound¬ 
ing  was  found  to  be  capable  of  measuring  the 
thickness  of  sea  ice  from  thin  lead  ice  to  thick  mul¬ 
tiyear  ice.  As  a  result  of  the  1 985  study,  a  follow-on 
effort  to  improve  the  AEM  ice  measurement  system 


was  determined  to  be  desirable.  The  goals  of  the 
1986-87  effort  were  to 

1 .  Improve  the  processing  software  in  terms  of 
speed,  accuracy  and  ease  of  use. 

2.  Perform  system  optimization  studies. 

3.  Make  improvements  to  the  calibration  proce¬ 
dure. 

4.  Investigate  the  possibility  of  doing  real-time 
data  processing. 

5.  Design,  construct  and  test  a  3 .5-m-long  an¬ 
tenna  (bird). 

6.  Design,  test  and  install  a  50-kHz  coil  in  the 
bird  for  the  purpose  of  improving  ice  thickness 
measurement  and  with  the  goal  of  measuring  ice 
conductivity. 

7.  Assemble  and  demonstrate  an  in-field  data 
processing  system. 

8.  Perform  a  field  validation  study  of  the  in¬ 
tegrated  AEM  hardware  and  processing  system. 

Advances  were  made  in  most  of  the  above 
areas,  particrolarly  in  improved  in-field  processing 
and  the  development  of  a  small  AEM  bird.  In  add¬ 
ition,  funds  were  provided  to  the  University  of 
CalifomiaatBerkeleyfordevelopmentofimproved 
AEM  interpretation  techniques  related  to  sea  ice 
thickness  and  conductivity  (Becker  et  al.  1987). 


PRE-FIELD-DEPLOYMENT  ACTIVITIES 

Software  development 

The  software  developed  for  processing  and 
analyzing  AEM  data  for  the  determination  of  sea 
ice  thickness,  seawater  conductivity  and  depth 
(Kovacs  et  al.  1987a)  was  improved  to  permit  fast¬ 
er  turnaround  of  the  raw  data  into  interpreted  pro¬ 
file  results.  The  processing  software  package  used 
to  analyze  the  1985  data  consisted  of  about  10  pro¬ 
grams,  whereas  the  current  one  used  four:  a  tape- 
to-disk  reading  program,  a  "master  calibration" 
program  called  MASCAL,  an  inversion  routine 


called  PRICE87,  and  a  plotting  routine  called 
LINPLOT.  The  entire  process  of  reading  the  raw 
AEM  data,  correcting  for  drift,  scaling  of  the  data, 
mterpreting  the  reduced  data  and  plotting  the 
results  has  been  reduced  to  a  few  hours  per  flight. 
Tl\is  compares  to  one  or  two  days  per  flight  with 
the  software  previously  used  to  analyze  similar 
data  collected  in  1985.  Data  processing  is  still  the 
slowest  single  operation,  with  a  throughput  of 
roughly  one  data  point  per  second.  This  time  in¬ 
cludes  determination  of  ice  thickness,  ice  conduc¬ 
tivity  and  water  depth.  Determination  of  these 
parameters  plus  water  conductivity  and  bottom 
conductivity  can  require  perhaps  10  seconds  per 
sounding  point.  When  only  ice  thickness  informa¬ 
tion  is  desired,  interpretation  could  proceed  using 
only  the  high  frequency  data.  This  can  speed  up 
the  interpretation  process  by  more  than  a  factor  of 
three. 

System  optimization 

Computer  studies  for  a  3.5-m-long  bird  indi¬ 
cated  that,  if  the  antenna  system's  noise  level  were 
on  the  order  of  0.1  ppm  of  the  primary  magnetic 
field  strength  at  the  receiver  (comparable  to  1  ppm 
achievable  in  the  larger  7.5-m-long  bird  used  in 
1985),  the  small  bird  would  be  able  to  resolve  ice 
cond  uctivity  down  to  about  0.01  S/  m.  Tlie  50-kHz 
quadrature  data  are  the  most  applicable  to  this  ice 
parameter  determination.  However,  theestimation 
of  ice  conductivity  is  very  sensitive  to  systematic 
errors  such  as  drift,  improper  calibration,  and 
laser  altimeter  errors,  especially  at  frequencies  of 
less  than  100  kHz.  A  working  frequency  of 200-300 
kHz  is  now  considered  necessary  for  good  ice  con¬ 
ductivity  measurements  of  sea  ice  (Becker  et  al. 
1987),  which  exhibits  a  wide  conductivity  range 
dependent  upon  the  ice  thickness,  salinity,  tem¬ 
perature,  age,  etc.  (Kovacs  et  al.  1987b). 

An  optimal  frequency  domain  ice-measurement 
system  could  utilize  three  coil  pairs  operating  at 
roughly  10  kHz,  50  kHz  and  300  kHz,  with  the  last 
two  coil  pairs  operating  in  the  vertical  coaxial 
mode  to  reduce  system  footprint  size  to  about  1.3 
times  the  bird  elevation  and  provide  very  good  sea 
ice  thickness  sounding  data.  The  system  footprint 
is  defined  as  the  area  at  the  ice/seawater  interface 
that  contributes  to  the  determination  of  the  average 
bird  height  above  the  sea  water.  The  effect  of  the 
footprint  is  to  smooth  out  the  seawater  depressions 
related  to  ice  draft  variation  and  in  turn  to  smooth 
out  the  ice  relief  with  what  amounts  to  a  moving 
average  filter  with  a  width  of  about  1.3  times  the 
bird  height  for  the  coaxial  coils  and  about  2  times 


the  bird  height  for  the  horizontal  coplanar  coils. 
The  1 0-kHz  coil  pair  should  give  good  bathymetry 
to  perhaps  2  m  below  the  ice,  which  would  be  suffi¬ 
cient  to  identify  grounded  ice  and  shoals.  The  two 
higher  frequencies  would  provide  good  definition 
of  ice  thickness  and  ice  conductivity  respectively. 
An  overall  reduction  in  weight  to  about  100  kg, 
from  the  current  small  bird  weight  of  about  1 50  kg, 
is  also  within  reason  and  should  make  a  fixed- 
wing  aircraft  deployment  of  the  system  more  feas¬ 
ible  as  well. 

Improvements  to  the  calibration  procedure 

System  calibration  proceeded  c  n  three  levels. 
On  the  basic  level,  an  investigation  of  the  accuracy 
of  the  standard  (^-coil  and  ferrite  bar  calibration 
method,  as  used  by  mineral  prospecting  firms, 
yielded  the  result  that  the  phase  calibration  is 
more  accurate  (even  at  high  frequencies)  than  the 
amplitude  cahbration.  However,  the  amplitude 
calibrations  were  not  in  agreement  with  the  calcu¬ 
lated  ones  determined  at  a  known  elevation  over  a 
surface  or  mass  of  known  conductivity  — or,  in 
other  words,  with  known  ground  truth.  An  effect 
of  the  secondary  field  produced  by  the  ground 
during  ground  calibration  probably  contributed 
to  the  disagreement  between  the  calculated  and 
measured  amplitude  calibrations.  It  had  been 
thought  that  the  small  but  nonzero  conductivity  of 
the  ferrite  bar  used  for  phasing  was  a  limiting 
factor,  but  this  effect  did  not  appear  at  frequencies 
up  to  50  kHz.  In  contrast,  the  standard"external 
Q-<oir'  used  for  ground  calibration  of  the  EM  sys¬ 
tem  appears  to  run  into  serious  difficulty  at  high 
frequencies,  mainly  due  to  the  parasitic  capacitance 
of  the  coil.  A  small-diameter  coil  was  constructed 
which  possessed  a  much  higher  self-resonant  fre¬ 
quency,  and  gave  much  better  results  than  did  the 
"standard"  (^-coil. 

The  second  level  of  investigation  centered  on 
the  design  and  construction  of  an  "absolute"  cali¬ 
bration  technique  that  could  be  used  in-flight.  This 
would  eliminate  the  secondary  field  effects  asso¬ 
ciated  with  on-the-ground  system  calibration.  A 
novel  technique  based  on  increasing  the  area  of  the 
receiver  coil  in  a  precise  way  was  devised  and 
tested.  However,  it  did  not  operate  as  required .  In¬ 
teraction  between  the  calibration  coil  and  the  re¬ 
ceiver  coil  caused  severe  departures  of  the  calibra¬ 
tion  signal  from  correct  levels,  especially  at  high 
frequencies.  As  it  stands,  this  technique  is  not  a 
satisfactory  stand-alone  calibration  aid. 

The  final  level  of  calibration  investigation  in¬ 
volved  the  automation  of  the  "ground-truth  re- 
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calibration"  process.  In  this  process,  the  AEM  bird 
altitude  and  laser  altimeter  results  obtained  over  a 
site  with  known  ice  thickness,  water  depth,  and 
water  conductivities  were  compared  with  the  com¬ 
puter-generated  responses  for  the  same  bird  posi¬ 
tion  and  orientation  over  the  site.  From  these  com¬ 
parisons,  corrections  to  the  scaling  factors  for  the 
AEM  data  were  derived.  These  scaling  factors  re¬ 
late  the  drift-corrected  EM  data,  in  bits,  to  the  final 
EM  results  in  parts  per  million  of  the  primary  field 
at  the  receiver.  Tlie  raw  data  were  then  recalibrated 
with  the  revised  scaling  factors  for  reinterpretation 
to  bring  the  AEM  results  in  line  with  the  measured 
ground  truth  data. 

It  should  be  emphasized  that  the  system's  cali¬ 
bration  is  a  static  quantity,  and  that  once  a  good 
calibration,  based  on  ground  truth,  has  been  estab¬ 
lished  it  should  be  stable  for  long  periods  of  time 
(months  or  even  years),  as  long  as  the  coils  and 
electronics  of  the  bird  and  the  receiver  console  par¬ 
ameters  are  not  altered.  Thus,  the  first  effort  on  a 
survey,  after  adjustments  have  been  made  to  the 
bird  or  console,  would  be  to  establish  a  good  cali¬ 
bration  from  ground  truth.  When  an  absolute  cali¬ 
bration  technique  has  been  fully  developed  and 
implemented,  this  step  will  not  be  required. 

Real-time  processing 

Substantial  effort  went  into  investigation  of  the 
feasibility  of  real-time  interpretation  of  AEM  ice 
measurement  data.  Two  factors  have  stymied 
efforts  to  process  AEM  data  in  real  time  to  yield  sea 
ice  thickness  or  other  desirable  quantities.  These 
are  the  effects  of  system  drift  and  the  difficulty  of 
performing  the  complex  computations  required  in 
a  small  reliable  package  suitable  for  airborne  oper¬ 
ations.  The  software,  which  was  configured  to  run 
on  a  ground-based  MicroVAX  II  computer,  can 
handle  from  1-3  samples  per  second.  However, 
improvements  in  the  efficiency  of  existing  software, 
as  well  as  mathematical  developments  by  several 
research  groups  in  the  past  year,  offer  the  possibility 
of  interpretation  rates  of  up  to  50  points  per  second 
or  about  one  data  point  per  meter  at  180  knots. 
Even  if  these  approaches  prove  unsatisfactory, 
computer  technology  has  progressed  so  rapidly 
that,  with  a  new  digital  receiver  incorporating  a 
80386-based  computer  coupled  with  a  small  array 
processor,  it  should  be  possible  to  process  AEM 
data  in  real  time  at  the  data  acquisition  rate. 

The  key  factor  limiting  progress  in  real-time 
processing  is  thus  not  the  interpretation  system 
but  the  drift  characteristics  of  present  AEM  systems. 
In  particular  the  problem  is  the  presence  of  so- 


called  "nonlinear"  drift:  instrumental  drift  that 
departs  significantly  from  a  linear  time  depen¬ 
dence.  Drift,  for  the  purposes  of  this  report,  can  be 
defined  as  the  change  vs  time  that  would  be  ob¬ 
served  in  the  AEM  system's  outputs  at  high  alti¬ 
tude,  when  there  is  no  change  occurring  in  the  ac¬ 
tual  electromagnetic  secondary  field  because  there 
are  no  conductive  structures  nearby  to  host  eddy 
current  flow.  If  this  drift  were  highly  linear  as  a 
function  of  time,  it  would  be  easily  removed  by 
making  high-altitude  measurements  at  both  ends 
of  a  survey  flight  line  and  then  interpolating  be¬ 
tween  these  "zeroing"  points.  A  combination  of 
good  hardware  design,  assembly  and  testing 
should  have  virtually  eliminated  nonlinear  drift. 
Unfortunately,  these  aspects  were  not  adequately 
addressed  by  the  manufacturer  but  will  be  in¬ 
corporated  into  the  next  phase  of  our  AEM  sea  ice 
measurement  program. 

Small  bird  design  and  implementation 

A  prototype  AEM  antenna  bird,  incorporating 
limited  bathymetric  capability  in  addition  to  its 
primary  design  goal  of  sea  ice  thickness  measure¬ 
ment,  was  constructed  by  Geotech  Ltd.  (Fig.  1). 
The  bird  was  suspended  by  a  Kevlar-cored 
"strength  member"  from  the  survey  helicopter. 
The  principal  bird  characteristics  were  a  coil  sep- 


Figurel.  Standard  size  AEM  bird  (a)  used  in  1985  and 
new  shorter  bird  (b)  used  in  the  1987  program.  The 
relative  size  difference  can  be  observed  bp  comparison 
with  the  length  of  the  spreader  bar  with  drag  ski  rt  which 
is  the  same  in  both  photographs. 


3 


aration  of  3  m  and  a  combined  bird  and  tow  cable 
weight  of  120  kg,  or  less  than  half  the  weight  of  a 
conventional  7-  to  8-m-long  bird  and  tow  cable 
(about 255  kg).  Since  the  smaller  coil  separation  re¬ 
duces  the  amplitude  of  the  secondary-field  anom¬ 
aly  from  a  given  conductor  at  a  given  standoff  dis¬ 
tance,  it  was  necessary  to  minimize  all  sources  of 
EM  noise  produced  within  the  bird.  Bird  EM  noise 
arises  from  a  variety  of  sources:  movement  of  the 
coils  with  respect  to  one  another  (due  to  bird  bend¬ 
ing  and  thus  internal  coil  platform  flexure  or  too 
loose  fasteners),  movement  of  conductors  (such  as 
bird  wiring)  and  movement  of  electronic  housings 
and  heat  sink  components  within  the  bird.  Sferic 
and  cultural  electromagnetic  no^e  picked  up  in 
the  coils  or  bird  wiring  also  contributes  to  the  over¬ 
all  measured  no'se  level.  To  achieve  a  signal-to- 
noise  ratio  in  a  small  bird  with  a  3-m  coil  separation 
that  is  comparable  to  that  of  a  conventional  AEM 
bird  with  a  6-  to  7-m  coil  separation,  the  noise  level 
must  be  about  1 0  times  lower.  This  noise  reduction 
was  addressed  by  making  the  bird  shell  more  rigid 
through  the  use  of  Kevlar-epoxy  materials,  and 
using  substantially  larger  dipole  moments  in  the 
transmitters  than  is  usual.  Lower  noise  preampli¬ 
fiers  were  also  fabricated. 

A  3-m  coil  spacing  was  chosen  for  the  approx¬ 
imately  3.5-m-long  bird  with  two  goals  in  mind: 
first,  to  minimize  the  difficulty  of  shipping  the  bird 
to  remote  locations  (a  major  expense  for  birds  7  to 
8  m  long),  and  second,  to  take  the  first  step  toward 
making  a  bird  that  could  be  towed  by  a  fixed-wing 
aircraft. 

The  coil  configurations  used  in  the  small  bird 
were  horizontal  coplanar  for  the  811 -Hz  (FI)  and 
49,927-Hz  (F3)  frequency  coil  pairs  and  coaxial  for 
the  4511 -Hz  (F2)  frequency  coil  pair.  The  coaxial 
coil  arrangement  was  selected  to  take  advantage 
of  the  inherently  narrower  "footprint"  of  this  con¬ 
figuration  and  thereby  obtain  data  for  enhancing 
interpretability  of  "two-dimensional"  features  like 
pressure  ridges. 

The  performance  of  the  small  bird  on  the  ground 
and  in  the  air  was  encouraging.  The  noise  level  for 
the  quietest  coil  pair  (4.5  Hz)  using  a  0.1-second 
detector  time  constant  was  approximately  0. 1  ppm, 
which  is  comparable  to  the  1-ppm  noise  level  of  a 
conventional  AEM  bird  with  a  7-m  coil  separation. 
Drift  was  on  the  order  of  5  ppm/hour  for  the  same 
coil  pair,  which  is  worse  than  expected  although 
quite  manageable.  Noise  and  drift  levels  were 
higher  in  the  other  channels,  particularly  the  50 
kHz.  The  resulting  errors  affected  the  quality  of 
the  field  data.  Thorough  test  and  evaluation  of  the 


bird  and  receiver  electronics  by  the  manufacturer, 
prior  to  slripment  to  Alaska,  should  have  uncovered 
this  noise  problem.  The  major  noise  sources  were 
found,  but  only  after  completion  of  the  Prudhoe 
Bay  field  test  program.  Drift  in  all  data  channels 
was  primarily  a  problem  caused  by  poor  stability 
in  the  analog  receiver  console.  Again,  this  problem 
should  be  alleviated  by  use  of  a  digital  receiver  or, 
at  the  very  least,  by  use  of  a  properly  stabilized  an¬ 
alog  console.  Replacement  of  the  transmitter  mod¬ 
ules  with  current-regulated  transmitters  should 
further  reduce  drift  in  the  bird's  electronic  system. 
With  the  above  modifications,  new  equipment, 
etc.,  further  shortening  of  the  bird  through  reduc¬ 
tion  of  the  coil  separation  may  be  feasible. 

Other  improvements 

A  commercial  16-bit  digital  data  recording  sys¬ 
tem  was  used  for  acquisition  of  all  survey  data. 
This  system  was  an  improvement  over  the  rather 
primitive  1 2-bit  data  acquisition  system  used  dur¬ 
ing  the  1985  survey .  These  recording  systems  gather 
analog  and  digital  inputs  from  a  variety  of  sources 
and  store  them  with  timing  information  on  mag¬ 
netic  tape.  The  recording  system  did  develop  hard¬ 
ware  problems  during  the  field  survey  that  caused 
a  two-day  delay  and  much  concern.  A  field  repair 
was  successful. 

A  video  flight  path  recovery  system  was  used 
that  gave  good  picture  quality  of  the  terrain  over¬ 
flown  and  allowed  for  digital  annotation  of  the 
video  record.  This  system  incorporated  a  camera, 
with  a  wide-angle  lens  and  image  sensor,  a  video 
annotation  board  and  a  compact  video  tape  re¬ 
corder.  The  use  of  a  microphone  for  the  operator  to 
record  cockpit  conversation  on  the  audio  track  of 
the  video  tape  cartridge  made  analysis  of  the  video 
record  easier. 

A  global  positioning  system  (GPS)  was  used  for 
in-flight  navigation  when  possible;  useful  three- 
satellite  "visibility  windows"  lasted  for  about  5 
hours /  day  in  1 987.  Positions  in  latitude  and  longi¬ 
tude  were  recorded  by  the  data  system  when 
available  and  served  to  guide  reconstruction  of 
flight  paths.  The  GPS  antenna  was  mounted  on  the 
helicopter  roof  above  the  cockpit.  This  position 
was  not  expected  to  work  very  well,  but  was  in  fact 
satisfactory;  the  only  problem  encountered  was 
occasional  blocking  of  the  signal  from  satellites 
low  on  the  horizon  by  the  engine  housing. 

A  new  laser  altimeter  was  installed  in  the  bird  to 
replace  the  less  accurate  unit  used  during  the  1985 
field  trials.  The  new  altimeter  is  much  more  accu¬ 
rate  (±1  cm)  than  the  old  one  (±20  cm).  Difficulties 
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Figures.  First-year  pres¬ 
sure  ridge  site.  The  arrows 
indicate  the  ends  of  the 
230-m-long  grid  later  es¬ 
tablished  at  this  location. 
Structure  in  background 
is  a  bottom  founded  con¬ 
crete  drilling  platform. 


Figure  2.  Outline  of 
50-m-wide  by  450-m- 
long  grid  area  on  mul¬ 
tiyear  floe  1.  Arrows 
point  to  grid  center- 
line.  Foot  traffic  during 
the  drilling  operation 
caused  contrast  seen  in 
the  snow  surface. 


were  encountered  in  programming  the  control 
console  for  this  altimeter.  This  appeared  to  be 
linked  to  problems  encoimtered  in  altimeter  data 
timing,  as  will  be  discussed  later. 


FIELD  ACTIVITIES 

Snow,  ice  and  seawater  data 

Before  making  the  AEM  sounding  flights,  four 
sites  were  selected  within  the  fast  ice  area  near 
Prudhoe  Bay,  Alaska.  At  each  site  extensive 


"ground  truth"  information  was  collected.  Tlrree 
sites  were  on  second-year  ice  floes,  hereafter  called 
multiyear  ice,  and  one  site  comprised  a  first-year 
pressure  ridge  surrounded  by  first-year  sea  ice. 

On  each  multiyear  ice  floe  (Fig.  2)  and  at  the 
first-year  pressure  ridge  (Fig.  3),  a  grid  was  laid 
out  with  stations  spaced  5  m  apart.  Grid  width  was 
50  m.  Grid  lengths  varied  and  were  450,  200,  200 
and  230  m  for  multiyear  floes  1,  2  and  3  and  the 
first-year  ridge  site,  respectively.  At  each  grid  sta¬ 
tion,  the  snow  and  ice  thicknesses  were  detennined, 
as  were  the  ice  freeboard  and  keel  depth.  Two  oil- 
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fired  hot  water  drilling  systems  were  used  to  drill 
most  of  the  holes.  About  75  holes  were  drilled 
using  a  hand-held  electrically  driven  6-cm-diam 
continuous  flight  auger  system.  In  all,  about  2,500 
holes  were  drilled. 

The  purpose  of  the  grids,  and  related  drill  hole 
measurements,  was  to  define  the  ice  thickness  dis¬ 
tribution  at  each  site.  The  data  were  used  to  evaluate 
if  the  snow-ice  thickness  determined  by  the  AEM 
system  agreed  with  the  measured  values.  In  our 
previous  work  (Kovacs  et  cil.  1 987a),  we  determined 
that  the  diameter  of  the  AEM  footprint,  or  the  area 
over  which  the  integrated  snow-ice  thickness  is 
estimated,  ranges  from  one  to  two  times  the  bird 
altitude.  Therefore,  if  the  bird  is  flown  25  m  above 
the  sea  ice,  then  the  footprint  is,  as  a  first  approxi¬ 
mation,  at  least  25  m  in  diameter. 

The  object  of  the  AEM  ice  thickness  study  was 
to  fly  down  the  center  of  each  grid,  sounding  ice 
thickness  and  referring  these  measurements  to  the 
location  of  markers  placed  at  25-m  intervals  on  the 
ice  surface.  Common  fiducial  numbers  recorded 
on  the  AEM  data  and  the  flight  path  video  allowed 
for  cross  correlations  of  the  ice  marker  location  vs 
the  AEM  data.  Past  experience  showed  that  it 
would  be  extremely  difficult  to  fly  the  bird  down 
a  narrow  track  laid  out  on  the  ice.  Some  bird  wan¬ 
dering  was  expected  due  to  helicopter  movement. 
To  compensate  for  this,  the  grid  was  made  wider 
than  the  anticipated  EM  footprint. 

An  additional  reason  for  the  wide  grid  was  to 
assess  the  track  width  needed  to  provide  average 
snow-ice  thickness  information  that  is  representa¬ 


tive  of  the  integrated  footprint  area  thickness  de¬ 
termined  by  AEM  sounding. 

In  addition  to  the  four  grid  sites,  three  long 
sounding  lines  were  established  across  the  fast  ice 
with  the  use  of  the  satellite  GPS.  At  random  dis¬ 
tances  along  each  line,  ice  stations  were  set  up  at 
which  snow  depth,  ice  thickness  and  water  depth 
were  measured.  At  several  stations,  seawater  tem¬ 
perature  and  conductivity  were  measured  and  ice 
cores  taken  for  ice  property  determinations.  The 
above  lines  were  later  flown  with  the  AEM  system 
with  the  goal  of  correlating  the  AEM  measured  re¬ 
sults  and  the  station  ground  truth  measurements. 

Another  site  selected  for  verification  purposes 
was  a  lead  containing  ice  0.18  to  0.22  m  thick.  The 
new  ice  did  not  have  a  snow  cover. 

The  locations  of  flight  lines  containing  stations 
A,  B,  and  D  through  G  are  shown  in  Figure  4.  The 
flight  line  containing  stations  H  through  O  is  shown 
in  Figure  5,  as  is  the  general  location  of  multiyear 
floe  study  sites  1, 2,  and  3.  The  first-year  pressure 
ridge  site  was  located  to  the  northwest  of  Reindeer 
Island  and  the  refrozen  lead  site  was  northeast  of 
this  island.  Right  line  station  positions  were  also 
determined  by  use  of  the  satellite  GPS.  Table  1 
gives  station  latitude  and  longitude  as  well  as  each 
site's  average  snow  and  ice  thickness,  and  the 
temperatures,  salinity  and  depth  of  water  under 
the  ice. 

Ice  thickness,  over  which  AEM  sounding  flights 
were  made,  varied  from  a  few  centimeters  to  about 
21  m  in  ridges.  Under-ice  water  depths  ranged 
from  zero  at  grounded  fast  ice  sites  to  over  30  m. 


Table  1.  Station  GPS  location  and  related  snow,  ice  and  water 
ground  truth  data. 


Thkkiiesii  Wahr 


Station 

Lititmk 

Longitude 

Sttow 

(m) 

Ice 

(m) 

Depth* 

<m) 

salinity 

(%r)' 

Temp. 

(Xi) 

A 

70‘’20.2'N 

148''20.4'W 

0.28 

1.70 

0.55 

70 

3.8 

B 

70°22.r 

148°21.6' 

0.25 

1.19 

Grounded  ice 

D 

70°24.3' 

148°28.8' 

0.16 

1.77 

1.95 

31.6 

1.75 

E 

70‘’24.7' 

148''28.0' 

0.14 

1.75 

4.4 

31.8 

1.75 

F 

70°26.r 

148°24.0' 

0.30 

1.68 

5.05 

30.9 

1.7 

G 

70°27.r 

148°22.8' 

0.17 

1.82 

5.1 

31.5 

1.75 

H 

70°30.4' 

149’>12.4' 

0.10 

1.80 

0.45 

40 

2.5 

I 

70<’32.4' 

I49W.0' 

0.15 

1.76 

9.1 

— 

— 

1 

70°33.0' 

149'’04.8' 

0.15 

1.73 

10.0 

31.0 

1.7 

K 

70°33.4' 

149”04.5' 

0.10 

1.74 

8.85 

31.3 

1.74 

L 

70°34.4' 

149°02.2' 

0.15 

1.78 

12.1 

— 

M 

70°35.r 

148'’59.7' 

0.15 

1.80 

15.15 

— 

— 

N 

70°36.7' 

148°55.6' 

0.05 

0.98 

14.1 

— 

— 

O 

70°37.4' 

148°53.0' 

0.05 

0.89 

18.7 

— 

— 

*  Depth  under  ice. 
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Figure4.  Location  of  flight  lines  F15L1  and  F12L1  and  ground  truth  stations  A, Band  D-G. 


Water  temperature  and  conductivity,  as  deter¬ 
mined  by  instrumentation  lowered  below  the  ice 
sheet,  were  typically -1 .75  C  and  2.6  S/m,  respec¬ 
tively,  at  locations  outside  of  confined  bays  and 
lagoons.  In  Prudhoe  Bay  at  station  A  (Fig.  4),  high¬ 
er  salinity  water  with  a  conductivity  of  about  4.6 
S/m  occurred,  and  in  Simpson  Lagoon  south  of 
Cottle  Island  at  station  H  (Fig.  5)  the  water  conduc¬ 
tivity  was  about  3.4  S/m.  The  higher  conductivity 
of  the  seawater  at  stations  A  and  H  occurs  because 
these  bodies  of  water  become  isolated  after  the  sea 


ice  thickens  and  becomes  grounded  on  the  sur¬ 
rounding  shallows.  No  longer  flushed  away  by 
currents,  the  brine  expelled  into  the  seawater  dur¬ 
ing  ice  growth  increases  the  salinity  and  therefore 
the  conductivity  of  the  confined  water.  This  en¬ 
richment  also  depresses  the  freezing  point  of  the 
water,  as  indicated  in  Table  1. 

AEM  surveys 

EKxring  the  course  of  the  field  experiment,  nearly 
40  passes  were  made  over  the  survey  areas  previ- 
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Figiirc5.LocatkviofflightliueF12L2,slatiousH  through  Oaud  the  three luulthjcaricc floes 
at  which  AEM  soundings  ivcre  made. 


ously  discussed.  The  final  flight  consisted  of  a 
northing,  a  long  flight  out  over  the  pack  ice.  Repeat 
passes  were  made  over  all  sites  for  which  ground- 
truth  information  had  been  obtained.  Most  flying 
was  done  over  the  four  grid  sites.  Total  line  coverage 
was  on  the  order  of  100  km. 

Every  attempt  was  made  to  limit  survey  altitude 
in  an  effort  to  reduce  footprint  size.  At  the  four 
grid  sites,  multiple  passes  were  made  to  ensure 
that  three  or  more  "good  runs"  down  the  center  of 
the  grid  would  be  obtained.  Rights  were  made  un¬ 
der  GPS  control,  when  three  or  more  "visible"  sat¬ 
ellites  made  this  possible,  along  the  survey  lines 
with  ground  truth  stations,  the  grid  sites,  and  the 
northing  line  and  at  several  other  locations  .  The 
satellite-positioned  flight  lines  discussed  in  this 
report  are  shown  in  Figure  6.  Tire  bird  was  also 


taken  to  zeroing  altitude  (about  3(X)-m  elevation) 
before  and  after  each  sounding  run.  At  this  altitude 
the  energy  emitted  by  the  transmitter  coils  is  no 
longer  capable  of  inducing  an  eddy  current  sec¬ 
ondary  system  in  the  conductive  seawater  that  is 
detectable  at  the  receiver  coils.  Thus,  no  secondary 
electromagnetic  field  is  sensed  by  the  receiver 
coils.  Measurement  of  baseline  response  of  the 
AEM  system  before  and  after  each  sounding  run 
helped  to  determine  the  amount  of  drift  and  al¬ 
lowed  for  the  removal  of  this  drift  from  the  data.  In 
almost  all  cases,  this  approach  to  baselining  paid 
off  in  terms  of  simplified  post-flight  data  processing 
and  improved  data  quality.  However,  in  some  cas¬ 
es,  nonlinear  drift  could  not  be  compensated  for, 
particularly  on  long  lines. 


Fii^iire  6.  Location  map  showing 
flight  lines  determined  bp  use  of  the 
satellite  global  positioning  si/stem 
‘  (CPS). 


Calibration 

System  calibration,  initially  performed  using 
standard  external  on-ground  techniques,  was  re¬ 
fined  by  flying  over  a  known  thickness  of  sea  ice. 
Tlie  standard  groimd  calibration  procedure,  under¬ 
taken  after  system  warmup,  consisted  of  first  align¬ 
ing  the  phase  of  the  system  ("phasing"),  then  set- 
tingitsgain.  Phasing  was  accomplished  by  rotating 
a  ferrite  bar  in  and  out  of  alignment  with  the  max¬ 
imum  electromagnetic  field  of  a  given  coil  pair. 
Tliis  rotation  genera  ted  a  purely  in-phase  variation 


in  the  measured  secondary  field.  Any  deflection 
observed  in  the  out-of-phase  or  quadrature  signal 
amplitude  during  this  procedure  is  evidence  of  an 
incorrect  phase  adjustment  and  waseliminated  by 
successive  phasing  adjustments  at  the  receiver 
module  for  the  particular  coil  and  fret.]uency  being 
calibrated.  Tlie  gain  was  set  by  placing  an  "external 
Q-<oir'  on  a  test  jig  (Fig.  7)  at  a  known  distance 
from  the  receiver  coil  in  question.  A  contact  was 
then  opened  and  closed,  to  make  and  break  a 


Figure  7.  Ground  calibration  using  the  standard  Q-coil.  Black  line  above  bird  is  the 
runiimj  pavement. 
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hined  circuit  consisting  of  the  coil  and  a  calibration 
resistor.  This  circuit  generated  a  known  anomaly, 
with  equal  in-phase  and  quadrature  components, 
which  was  used  to  adjust  the  amplitude  response 
of  the  system. 

The  other  form  of  calibration  used  was  ground- 
truth  recalibration.  This  technique  is  far  more  accu¬ 
rate  than  the  standard  on-ground  calibration  tech¬ 
nique  when  good,  unambiguous  ground-truth  and 
altimeter  data  are  available.  However,  when  ice  of 
irregular  thickness  is  present  and/or  the  laser 
altimeter  is  not  performing  properly,  systematic 
errors  can  arise  from  this  procedure  that  will  bias 
the  results,  particularly  in  the  ice  thickness  estimate. 
The  ideal  situation  for  ground-truth  recahbration 
(assuming  proper  altimeter  operation)  is  therefore 
a  wide  area  of  deep  ice-free  water .  Since  such  lead 
conditions  may  not  be  available  when  required,  an 
internal,  absolute  calibration  system  will  be  essen¬ 
tial  for  a  fully  operational  AEM  sea  ice  measure¬ 
ment  system. 

AEM  data  processing 

A  field  processing  station  was  set  up.  It  consisted 
of  a  minicomputer  (MicroVAX  II),  two  terminals, 
a  cassette  tape  drive,  a  graphic  printer  and  related 
processing  software.  Analysis  of  the  AEM  data  fol¬ 
lowed  the  general  theory  outline  in  Kovacs  et  al. 
(1987a).  After  the  first  survey's  data  were  read  into 
the  computer  and  an  analysis  attempt  made,  serious 
difficulties  in  processing  began  to  be  encountered . 
This  was  characterized  by  an  inability  to  get  even 
moderately  reasonable  ice  thickness  results  out  of 
the  survey  data  while  at  the  field  station.  Resolution 
of  this  problem  is  discussed  below. 

POST-SURVEY  RESULTS 
Processing  problem 

The  source  of  the  problem  that  prevented  in¬ 
field  data  processing  eventually  turned  out  to  be 
surprisingly  simple:  a  time  offset  was  discovered 
between  the  laser  altimeter  data  and  the  rest  of  the 
survey  data.  The  altimeter  data  were  delayed  by 
about  2.5  seconds  from  the  AEM  data  on  the  stor¬ 
age  tapes.  This  was  a  serious  problem,  because 
normal  changes  in  laser  altitude,  due  to  flight 
height  and  bird  pitch  and  roll,  were  effectively  dif- 
ferentiated  and  added  into  the  ice  thickness, 
causing  very  substantial,  time-varying  errors  in 
the  interpreted  results.  Suffice  it  to  say  that  this 
time  skew  was  very  difficult  to  identify,  as  similar 
data  degradation  effects  might  have  been  caused 
by  a  variety  of  sources.  (Even  after  the  error  had 


been  found  to  be  in  a  time  skew  between  the  laser 
and  the  EM  measurements,  identification  of  the 
source  of  the  problem  proved  difficult.  Some  pos¬ 
sibilities  considered  were  errors  in  tape  writing  or 
transcription,  errors  in  MASCAL,  and  the  effects 
of  a  long  filter  time  constant  in  the  analog  receiver, 
as  described  by  Becker  and  Cheng  [1988].  Errors 
in  MASCAL  and  tape  transcripKon  were  eliminated 
as  the  cause  after  careful  analysis,  and  since  the 
time  constant  in  the  small  bird  analog  receivers 
was  0.1  second,  it  could  not  have  been  the  reason 
for  the  time  offset  either.)  In  the  final  analysis,  it 
appears  that  the  delay  arose  in  the  interface  between 
the  laser  process  and  the  data  acquisition  system. 
The  problem  can  be  eliminated  by  triggering  the 
laser  console  for  each  sample  desired. 

After  the  time  delay  was  accounted  for,  the 
AEM  data  were  recalibrated  using  the  snow  plus 
ice  (S-I)  thickness  measurements  from  the  flat  ice 
area  in  front  of  the  first-year  pressure  ridge  grid 
site.  It  should  be  noted  that  ground-truth  recali¬ 
bration  of  the  AEM  system  was  carried  out  at  this 
single  site,  and  that  this  calibration  was  then  used 
for  all  other  AEM  flight  data.  Profiles  of  the  S-I 
thickness  obtained  over  all  of  the  major  ice  struc¬ 
tures  were  then  prepared. 

Some  difficulties  were  encountered  during  the 
reduction  and  interpretation  of  the  AEM  data 
from  the  long  survey  lines.  Drift  in  the  F3  data  (50- 
kHz  frequency)  was  the  culprit  in  this  case.  Re¬ 
moval  of  the  F3  data  from  the  interpretation  solved 
this  problem. 

An  example  of  the  time-delay  corrected  data 
collected  during  flight  13,  line  3  (F13L3)  is  shown 
in  Table  2.  This  survey  line  consisted  of  a  pass 
down  the  200-m-long  grid  on  multiyear  floe  2.  The 
nine  reference  station  locations,  spaced  25  m  apart 
along  the  grid,  are  listed  by  fiducial  number  at  the 
bottom  of  the  table.  In  the  analysis  of  this  data  set 
the  conductivi  ties  of  the  sea  wa  ter  and  seabed  were 
set  at  2.6  and  0.5  S/m,  respectively,  rather  than 
being  interpreted  by  analysis  of  the  AEM  data. 
Note  that  the  bulk  ^I  conductivity  is  also  listed. 
These  S-I  conductivity  determinations  were  found 
to  be  inconsistent.  At  times  the  values  appeared 
very  representative  of  the  type  of  ice  over-flown 
(based  on  analyses  of  Kovacs  et  al.  1987b)  but  often 
they  varied  inconsistently  and  are  therefore 
considered  unreliable.  This  problem  appeared  to 
reside  with  the  nonlinear  drift  and  noise  experi¬ 
enced  with  the  AEM  system  and  not  with  the  algo- 
ritfim  used  to  estimate  the  bulk  S-I  conductivity.  In 
the  near  future  we  hope  to  rectify  this  problem 
through  hardware  improvements  and  the  use  of 
higher  frequencies. 


Table  2.  Example  of  processed  AEM  sounding  results.  Data  are  from  flight  line  F13L3. 
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Figure 8.  Suolu, freeboard audsmwplus ice thicknesscontourmapsof  multiyear floel  (courtesy  Anthony  Valentine). 


Ground  truth  data 

As  previously  stated,  large  grids  were  estab¬ 
lished  at  four  sea  ice  sites.  Their  50-m  width  was 
divided  into  11  longitudinal  lines,  A  through  K, 
spaced  5  m  apart.  Grid  length  was  similarly  di¬ 
vided  by  lateral  rows  5  m  apart.  At  the  intersections 
of  each  row  and  line,  the  snow  and  ice  thickness 
and  freeboard  were  determined  by  drill  hole  meas¬ 
urement.  These  data  were  used  to  construct  an  S-I 
thickness  contour  map  for  multiyear  floe  grid  sites 
1,  2  and  3.  These  maps  are  shown  in  Figures  8, 9 
and  10,  respectively.  The  relief  in  Figure  8  shows 
that  the  ice  was  thicker  at  each  end  of  the  grid,  as 
well  as  at  about  the  260-m  distance  location,  and 
that  ice  relief  features  tended  to  be  perpendicular 
to  the  long  axis  of  the  grid.  A  melt  pond  is  centered 
at  about  the  junction  of  line  I  and  the  340-m  dis¬ 
tance  position. 

The  S-I  thickness  contours  for  grid  site  2  show 
thick  ice  (a  pressure  ridge)  extending  from  line  A, 
at  about  the  60-m  distance  location,  to  line  K,  at 
about  the  85-m  distance.  Another  thick  ice  feature 
(a  ridge)  runs  from  about  line  I,  at  a  distance  of 
about  1 10  m,  to  about  lineC  near  the  200-m  distance 
location. 

The  contour  map  for  grid  site  3  shows  a  major 
ice  thickness  increase  (a  ridge)  extending  across 
the  grid  begiiming  in  the  area  of  the  1 40-m-distance 
location. 

From  the  drill  hole  measurements  the  average 
row  snow  and  ice  thickness  and  ice  freeboard  and 
draft  were  determined  for  grid  lines  A-K,  C-I  and 
E-G  at  each  grid  site.  These  values  are  listed  in 
Tables  3  through  6  for  multiyear  floes  1  through  3 
and  the  first-year  pressure  ridge  site,  respectively. 
Note  that  in  Table  6  the  average  ice  freeboard  on 
row  45  at  the  first-vear  ridge  was  a  negative  value. 
Indeed,  the  ice  in  this  area  was  depressed  by  snow 
loading  and  by  ridge  subsidence.  Drilling  through 
the  ice  along  this  line  resulted  in  upward  movement 
of  seawater  into  the  snow. 

At  the  bottom  of  each  table  is  listed  the  average 
value  for  each  column  of  data.  The  deviation  in  the 
average  ice  thickness  for  lines  A-K,  C-I  and  E-G 
for  any  grid  site  does  not  exceed  4%.  This  suggests 
that  ice  thickness  measurements  needed  to  be  ob¬ 
tained  only  along  three  of  the  lines  in  order  to 
assess  the  average  grid  area  ice  thickness. 

Table  7  lists  the  mean  and  standard  deviation  of 
the  snow  and  the  ice  thickness,  the  ice  freeboard 
and  draft  and  the  combined  snow  and  ice  thickness 
for  lines  A-K,  C-I  and  E-G  for  each  grid  site.  The 
maximum  deviation  of  the  mean  values  for  the 
snow-ice  thickness  data,  the  thickness  that  would 
be  measured  by  the  AEM  system,  is  again  very  low 


(less  than  3%)  for  each  group  of  lines.  This  again 
suggests  that  data  collected  along  lines  C-I  may  be 
used  to  assess  the  mean  snow,  ice,  S-I  thickness 
etc.,  for  the  entire  related  grid  area.  For  the  multi¬ 
year  grid  sites,  the  standard  deviation  of  the  snow 
thickness  and  the  ice  freeboard  values  is  seen  to  be 
comparable  to  the  mean  thickness.  This  is  as  ex¬ 
pected  since  the  floes  had  considerable  undulating 
relief  and  significant  variation  in  snow  drift  ac¬ 
cumulation.  Such  topographic  variation  is  even 
more  pronounced  at  the  first-year  pressure  ridge 
site,  as  indicated  by  the  standard  deviations  shown 
in  Table  7. 

While  the  mean  S-I  thickness  values  for  lines 
A-K,  C-I  and  E-G  for  each  grid  site  are  quite  sim¬ 
ilar,  there  was  variation  between  them  on  a  row- 
by-row  basis  or  with  distance  down  the  grid  axis. 
This  is  graphically  shown  in  Figures  1 1  through  14 
for  grid  sites  1  through  3  and  the  first-year  ridge 
site,  respectively.  Our  assessment  of  the  S-I  thick¬ 
ness  data  indicated  that  line  group  C-I  was  more 
representative  of  the  average  S-I  thickness  varia¬ 
tion  along  the  grid  and  was,  therefore,  used  to 
evaluate  the  AEM-determined  measurements  of 
S-I  thickness. 

To  show  how  the  S-I  thickness  can  vary  along 
individual  grid  lines,  we  include  Figures  1 5  through 

17  for  multiyear  flows  1  through  3,  resp»ectively. 
The  difference  between  the  lines  in  each  figure  as 
well  as  between  the  related  average  S-I  thickness 
lines  given  in  Figures  11  and  13  indicate  that  the 
variation  can  be  appreciable.  This  is  not  unexpjected. 

Cross  sections  constructed  from  the  drill  hole 
data  forcombined  lines  A-K,  C-I  and  E-G  for  mul¬ 
tiyear  grid  sites  1  through  3  are  shown  in  Figures 

18  through  20,  respectively.  These  cross  sections 
again  give  a  visual  indication  of  the  floes'  snow 
thickness  variation  as  well  as  the  variation  in  their 
freeboard  and  keel  relief. 

To  give  an  appreciation  of  the  ice  and  snow 
thickness  variation  along  each  grid  line  at  the  first- 
year  ridge,  cross-section  plots  were  made.  These 
are  presented  in  Figure  21  for  lines  A  to  K.  In  the 
area  of  the  grid,  the  keel  of  the  ridge  was  apparently 
in  contact  with  the  sea  bed.  Note  the  variations  in 
ridge  keel  and  sail  geometry  as  well  as  in  snow  ac¬ 
cumulation  depths. 

Averaging  the  cross  section  data  by  line  groups 
A-K,  C-I  and  E-G  gives  the  cross  sections  shown 
in  Figure  22.  The  ice  sheet  on  the  left  of  the  ridge 
had  a  relatively  uniform  snow-ice  thickness  of 
about  1.45  m.  The  ice  on  the  right  side  was  about 
15%  thicker.  Ice  blocks  in  the  ridge  varied  from  0.5 
to  1.5  m  in  thickness. 

It  should  be  noted  that  unlike  relatively  solid 
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freeboard  and  snow  plus  ice  thickness  contour  map  of  multiyear  floe  (courtesy  Anthony  Valentine). 


freeboard  and  snow  plus  ice  thickness  contour  map  of  multii/ear  floe  3  (courtesp  Anthonp  Valentine). 


Table  3.  Average  row  snow  depth  and  ice  thickness,  freeboard  and  draft  for  lines  A-K,  C-I  and  E-G  at  multiyear  floe  grid  site  1, 
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Table  4.  Average  row  snow  depth  and  ice  thickness,  freeboard  and  draft  for  lines  A-K, 
C-I  and  E-G  at  multiyear  floe  grid  site  2. 
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Table  5.  Average  row  snow  depth  and  ice  thickness,  freeboard  and  draft  for  lines  A-K,  C-I  and  E-G  at  multiyear  floe  grid  site  3. 
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Table  7.  Ice  and  snow  variation  in  the  mean  drill  hole  measured  data  between  line  groups  E-G, 
C-I  and  A-K  for  the  grids  at  the  three  multiyear  (MY)  ice  floes  and  the  first-year  pressure  ridge 
(FYR)  site. 


Grid 

site 

Grid 

lines 

Snow  thickness 

S-I  thickness 

Ice  thickness 

Ice-fri 

ebciard 

Ice  draft 

Mean 

(ml 

Std.  del', 
(ml 

Mean 

(m) 

Std.  del’, 
(ml 

Mean 

(ml 

Std.  del’, 
(ml 

Mean 

(ml 

Std.  del’, 
(ml 

Mean 

(ml 

Std.  dev. 
(ml 

MYl 

E-G 

0.115 

0.092 

3.945 

1.369 

3.830 

1.390 

0.380 

0.260 

3.450 

1.221 

MYl 

C-I 

0.122 

0.110 

3.934 

1.318 

3.812 

1.338 

0.391 

0.265 

3.421 

1.158 

MYl 

A-K 

0.125 

0.118 

3.958 

1.314 

3.833 

1.330 

0.397 

0.274 

3.436 

1.140 

MY  2 

E-G 

0.121 

0.144 

5.685 

1.341 

5.564 

1.378 

0.637 

0.515 

4.927 

1.025 

MY  2 

C-I 

0.124 

0.158 

5.639 

1394 

5.516 

1.434 

0.593 

0.470 

4.923 

1.106 

MY  2 

A-K 

0.111 

0.146 

5.505 

1.444 

5.394 

1.472 

0570 

0.456 

4.823 

1.155 

MY  3 

E-G 

0.077 

0.091 

4.067 

1.353 

3.989 

1.383 

0.437 

0.337 

3.553 

1.124 

MY  3 

C-I 

0.073 

0.094 

4.225 

1.506 

4.152 

1.536 

0.460 

0.341 

3.693 

1299 

MY  3 

A-K 

0.079 

0.109 

4.209 

1.543 

4.130 

1.578 

0.448 

0.334 

3.681 

1.341 

FYR 

E-G 

0.129 

0.141 

3.929 

4380 

3.801 

4.240 

0.412 

0.757 

3.389 

3.591 

FYR 

C-I 

0.120 

0.133 

3.935 

4.351 

3.816 

4.320 

0.433 

0.842 

3.382 

3.616 

FYR 

A-K 

0.111 

0.118 

3.955 

4.326 

3.844 

4.301 

0.412 

0.789 

3.432 

3.662 
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Figure  1 1 .  Average  snow  plus  ice  thickness  along  MY  floe  2 's  50-  x  450-ni  grid.  Profile  A-K  represents 
the  full  50-m  width  of  the  grid.  C-I  the  center  30-m  width  and  E-G  the  center  2  0-ni  width. 
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Figure  22.  Average  snoiv  plus  ice  thickness  along  MY  floe  2's  50-  x  200-nt  grid  for  szmth 
widths  A-K,  C-I  and  E-G. 
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Figure  13.  Average  snow  plus  ice  thickness  along  MY  floe  3's  50-  x  200-m  grid  for  swath 
widths  A-K,  C-I  and  E-G. 


Figure  14.  Average  snow  plus  ice  thickness  along  the  FY-ridge  50-  x225-m  grid  for  swath  widths  A-K,  C-I 
and  E-G. 


Figure  15.  Example  of  snow  plus  ice  thickness  variations  along  MY  floe  1  's  grid  for  lines  F  and  H  (see 
Fig.  6  for  line  location). 
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Figure  16.  Example  of  snow  plus  ice  thickness  variation  along  MY  floe  2' s  grid  for  lines  A  and 
K  (see  Fig.  7  for  line  location). 
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Figure!?.  Example  of  snow  plus  ice  thickness  variation  along  MY  floeS'sgrid  forlinesBand 
K  (see  Fig.  8  for  line  location). 


multiyear  pressure  ridges,  first-year  ridges  are  a 
conglomerate  of  ice  blocks  and  fragments.  Air  fills 
the  inter-block  voids  of  the  sail  and  seawater  fills 
the  voids  of  the  keel  extending  below  the  "solidifi¬ 
cation"  zone,  a  seawater  freezeback  zone  extending 
to  about  the  depth  of  the  adjoining  level  ice.  The 
EM  secondary  field  associated  with  this  high  con¬ 
ductivity  keel  structure  is  different  from  that  asso¬ 
ciated  with  the  low  seawater  or  brine  content  and 
therefore  lower  conductivity  of  a  multiyear  ridge 
keel.  Analysis  of  the  EM  data  taken  over  the  first- 
year  ridge  is  more  ambiguous  because  there  is  no 
well-defined  ice/water  interface  and  is  further 
complicated  by  the  relatively  narrow  width  and 
steep  slope  of  the  keel  (=  30  ),  its  three-dimensional 
geometric  variation  and  the  proximity  of  another 
conductive  zone,  the  sea  bed. 

Floe  freeboard  vs  thickness  analysis 

The  drill  hole  data  for  the  multiyear  floes,  listed 
in  Tables  3, 4  and  5,  were  used  to  assess  the  average 
density  of  the  ice  floes  and  their  equivalent  free¬ 


board  vs  draft.  The  equivalent  freeboard  is 
determined  by  converting  the  snow  cover  depth 
into  an  equivalent  ice  thickness  and  adding  this 
value  to  the  measured  ice  freeboard.  The  average 
snow  cover  was  of  unusually  low  density.  Indeed, 
almost  all  of  it  was  blown  away  during  a  period  of 
high  winds  at  the  end  of  our  field  program.  The 
average  measured  snow  specific  gravity  p^for  the 
snow  was  0.318  and  for  the  freeboard  ice  it  was 
0.880.  Therefore,  each  unit  of  snow  depth  is 
equivalent  to  0.364  unit  of  freeboard  ice.  For  all 
lines  (A-K)  on  multiyear  floe  1,  the  average  snow 
depth  was  0.125  m,  which  is  equivalent  to  0.046  m 
of  freeboard  ice.  The  average  measured  ice  free¬ 
board  was  0.397  m  and  the  equivalent  freeboard 
now  becomes  0.443  m.  The  average  ice  draft  D  was 
3.436  m,  which  gives  for  FJD  a  value  of  0.129  or  a 
ratio  of  1  to  7.756.  Buoyancy  considerations  dictate 
that 


5l+i  =  P«  (1) 

D  pi 
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Figure  21 .  Cross  section  along  FY  ridge  grid  for  lines  A~K.  Note  the  imriations 
in  the  ridge  keel  from  line  to  line  (5-m  spacing). 
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ion  along  FY  riiigc  grid  for  lines  A  through  K.  Note  the  variations  in  the  ridge  keel  from  hue 


Oltionc*  (ml 


where  p^,  is  the  specific  gravity  of  seawater  and  p. 
is  the  equivalent  specific  gravity  of  the  ice  floe.  The 
seawater  salinity  Sw  (in  parts  per  thousand)  was 
measured  as  31 .5%o.  Since  p  =  1  +  0.0008  Sw,  the 
specific  gravity  of  the  seawater  was  1 .025  and, 
and,  from  1,  pj  is  found  to  be  0.908.  The  apparent 
specific  gravity  of  the  submerged  ice  p^  may  be 
foimd  from 

Ps=Pw-Pi(VO>-  (2) 

Substituting  0.880  for  p^,  1 .025  for  p^^.and  0.129  for 
¥J D  in  eq  2,  a  value  of  0.912  for  p^  is  obtained. 

Making  similar  calculations  for  the  lines  A-K 
on  multiyear  floe  2  gives  0.1 26  for  F^/D  =  0.910  for 
p.  and  0.914  for  p  For  the  lines  on  multiyear  floe 
3  the  values  are  0.130  for  F^/D  ,  0.907  for  p^and 
0.911  for  p^.  The  variation  in  the  above  values  is 
due  in  part  to  ice  and  snow  property  variations  be¬ 
tween  the  ice  floes,  to  spatial  relief  variations  along 
each  grid  and  to  abrupt  Fy  D  variations  at  the  tran¬ 
sition  area  between  a  ridge  sail  and  the  surrounding 
ice  floe,  as  may  be  inferred  from  the  Kovacs  first- 
year  and  multiyear  pressure  ridge  models  (Wright 
et  al.  1981)  shown  in  Figure  23.  Large  pressure 


ridges  are  known  to  have  proportionally  less  dense 
ice  in  the  ridge  sails  and  increasing  ice  density 
with  depth  in  the  keel.  Such  density  variation  can 
upset  the  average  ratio  of  ice  freeboard  to  equivalent 
draft  of  -1:7.8  which  can  be  determined  from  the 
vs  D  values.  It  is  known,  for  example,  that  the 
mean  freeboard /draft  ratio,  obtained  from  sail 
height  and  keel  draft  measurements,  of  multiyear 
pressure  ridges  is  -1:3.2  as  shown  in  Figure  24. 
Ridges  in  the  floes  studied  may  have  biased  the 
equivalent  freeboard /ice  draft  results  in  favor  of 
the  lower  elevation  portion  of  the  ice  floes.  Tliis  is 
indicated  in  a  plot  of  the  average  grid  row  F^  vs  D 
values  which  are  listed  in  Table  8  and  plotted  in 
Figure  23.  The  data  in  Figure  25  indicate  that  for  ice 
drafts  over  about  6.5  m  there  is  a  marked  increase 
in  F^/D.  The  deep  draft  data,  as  presented  in  this 
figure,  are  simply  too  limited  to  support  this  view. 
The  analysis  for  the  line  passing  through  the  data 
in  Figure  25  did  not  include  the  four  data  points  for 
ice  drafts  over  6.5  m.  Because  of  this,  the  slope  of 
the  line  shown  (0.125)  is  lower  than  the  average  of 
0.128  for  the  three  F^/D.  values  previously  pre¬ 
sented.  Using  a  value  of  0.125  for  F  /D.and  p^,  of 
1.025  in  eq  1  gives  a  value  of  0.911  for  p^. 


Figure  23.  Kovacs  first-year  (a)  and  multiyear  (b)  pressure  ridge  models.  These  idealized  models  are  in 
hydrostatic  equilibrium  for  the  conditions  when  p=  0.92  and  p^  =  1 .025  for  seawater  of  salinity  =  31.5Vcc 
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Figure  24.  Multiyear  pressure  ridge  sail  height  vskeel  depth  (revised  from  Kovacs  1975  and 
1983). 
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Figure  25.  Multiyear  sea  ice  equivalent  freeboard  vs  draft. 


Tg.  Effective  Ice  Thickness  in?) 

Figure  26.  Multiyear  sea  ice  equivalent  freeboard  vs  effective  ice 
thickness. 


A  more  appropriate  presentation  of  the  data  is 
given  in  Figure  26.  Here  F^is  plotted  vs  the  effective 
ice  thickness  r_.  T  was  determined  by  adding  the 
equivalent  freeboard  ice  thickness  of  the  snow  to 
the  drill-hole-meas  xred  ice  thickness  from  Tables 
3, 4  and  5  (lines  A-K>  to  obtain  an  effective  total  ice 
thickness.  Therefore,  the  plot  in  Figure  26  should 
be  representative  of  snow-free  ice  floes.  The  linear 
regression  for  the  line  passing  through  the  data 
did  not  include  the  four  data  points  for  T,  values 
over  8  m  (for  reasons  previously  mentioned). 

The  slope  of  the  line  in  Figure  26  is  in  good 
agreement  with  a  regression  line  put  through  sim¬ 
ilar  F  ys  T^^data  by  Kovacs  (1977).  His  analysis  in¬ 
cluded  T  data  up  to  about  4  m,  as  shown  in  Figure 
27.  Beyond  this  thickness  thelimited  data  suggested 


an  abrupt  increase  in  the  F  ./T  , ratio.  The  data  pre¬ 
sented  in  Figure  26  indicate  that  the  abrupt  change 
shown  in  Figure  27  is  not  realistic.  It  may  be  noted 
that  the  data  in  Figure  27  fall  very  nicely  within  the 
data  presented  in  Figure  26. 

For  a  snow-free  multiyear  ice  floe,  the  analysis 
given  in  Figure  26  may  allow  for  a  reconstruction 
of  the  floe's  apparent  cross  section  from  the  AEM 
sounding  data.  Since  analysis  of  the  AEM  sounding 
data  gives  ice  thickness,  it  would  be  a  simple  mat¬ 
ter  to  use  the  equation  in  Figuiv  26  to  determine  an 
approximate  ice  freeboard. 

For  the  snow-covered  multiyear  sea  ice,  a  more 
representative  plot  of  the  measured  data  is  shown 
in  Figure  28.  This  plot  shows  the  measured  surface 
elevation  £  (snow  plus  freeboard)  vs  theS-I  thick- 
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Figure  27.  Multiyear  sea  ice  eijuivaleut  freeboard 
vs  effective  ice  thickness  (after  Ko'vacs  1977). 

ness  obtained  at  1896  drill  hole  sites  on  floes  1, 2 
and  3.  This  presentation  clearly  reveals  the  wide 
scatter  in  the  £  vs  S-I  data  due  to  the  snow  and  ice 
variations  previously  mentioned. 

The  linear  regressive  curve  shown  in  Figure  28 
wasdeterminedby  first  adding  increasingE  vsS-1 
data  until  a  linear  regression  curv'e  through  the 
data  intercepted  zero,  as  would  occur  in  nature. 
This  occurred  when  S-I  thickness  values  up  to  5.85 
m  were  included.  At  this  point  there  were  1593 
data  points  and  the  slope  of  the  line  was  0.1 24,  with 
a  correlation  coefficient  of  0.622  and  standard 
deviation  of  ±0.159.  Similar  results  were  obtained 
when  a  linear  regression  curve  was  forced  through 
zero  but  included  1 679  data  points  up  to  a  S-I  thick¬ 
ness  of  6.25  m.  This  linear  regression  curve  is  pre¬ 


sented  in  Figure  28  and  represents  an  E/S-I  thickness 
ratio  of  1  to  8.06,  which  from  eq  1  gives  a  value  of 
0.912  for  p.. 

The  third-degree  regression  curve  also  shown 
in  Figure  28  is  a  fit  to  the  entire  data  set.  An  appar¬ 
ent  cross  section  of  a  snow-covered  multiyear  ice 
floe  may  be  reconstructed  from  the  AEM-deter- 
mined  £  I  thicknesses  by  using  the  more  represen¬ 
tative  linear  curve  from  zero  up  to  a  S-I  thickness 
of  5.6  m  and  then  using  the  polynomial  equation 
for  higher  values  of  S-I  thickness.  However,  the  re¬ 
sulting  cross  section  may  be  of  limited  value  be¬ 
cause  of  the  low  correlation  coefficient  for  the 
curves  and  the  wide  scatter  of  the  data  shown  in 
Figure  28.  In  short,  the  resulting  cross  section 
would  not  reveal  the  highly  variable  hummocky 
relief  characteristic  of  multiyear  sea  ice  floes. 

AEM  sounding  results 

Station  K  (Fig.  5)  had  a  measured  S-I  thickness 
of  1 .79  to  1 .86  but  averaged  1 .84  m.  This  site's  S-I 
thickness  was  sounded  on  six  different  flights 
with  the  AEM  system.  The  bird  height  at  each 
sounding  was  between  19  and  23  m  above  the  sur¬ 
face.  On  average  the  soundings  were  made  at  21 .5 
m.  The  AEM  determined  S-I  thicknesses  were  1 .68, 
1.84,  1.90,  1.96,  2.45  and  2.50  m  which  gives  an 
average  of  2.06  m.  This  average  is  12%  higher  tlian 
the  average  measured  value.  On  an  individual  ba¬ 
sis  the  AEM  S-I  thickness  values  varied  from  9% 
below  to  36%  above  the  average  drill-hole-meas¬ 
ured  S-I  thicknesses.  Tliis  widespread  is  not  accept¬ 
able  and  must  be  improved  in  the  next-generation 
AEM  sea  ice  sounding  system.  The  spread  is  be¬ 
lieved  to  be  associated  with  nonlinear  system  drift. 

Nonlinear  drift  effects  may  be  seen  in  the  AEM 
S-I  thickness  results  (Table  9)  obtained  from  flight 


Figure  2S.  Siioiv  thickness  yins  ice  freeboard 
vs  total  snozv  pins  ice  thickness  for  ninitiyenr 
sea  ice.  (See  Appendix  A  for  bloivnp  of  this 
f/ginr.) 
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Table  9.  Drill-hole-measured 
vs  AEM-determined  snow- 
ice  thickness  for  stations 
H-O. 


Statioii 

Siiow-icc  thickm->i  (in) 

Drill  hole 

AEM 

H 

1.90 

1.15 

I 

1.91 

1.75 

J 

1.88 

2.37 

K 

1.84 

2.50 

L 

1.93 

2.24 

M 

1.95 

2.45 

N 

1.03 

0.79 

0 

0.94 

0.83 

lineF12L2  for  stations  H  and  O  (Fig.  5).  Table  9  lists 
the  average  drill-hole-measured  S-I  thickness  vs 
the  AEM  value  for  each  station  along  the  flight 
line. 

As  seen  in  Table  9  the  AEM  S-I  thickness  values 
at  the  beginning  of  the  flight  were  lower  than  the 
measured  values.  As  the  18-minute  flight  pro¬ 
gressed  the  AEM  values  gradually  increased  to 
where  they  exceeded  the  measured  ones  and  then 
decreased  again.  Time  between  zero  readings  was 
about  25  minutes.  Again  this  variation  is  unaccep¬ 
table  and  must  be  eliminated  before  AEM  sounding 
of  sea  ice  can  be  considered  a  viable  measurement 
technology. 

Due  to  the  short  flight  lines  and  therefore  the 
elapsed  time  between  the  zero  calibrations  made 
before  and  after  each  run,  drift  was  not  as  severe  a 
problem  in  the  AEM  data  collected  over  the  four 
grid  sites.  AEM  sounding  flights  typically  towed 
the  bird  at  an  elevation  of  23  ±3  m  above  the  sur¬ 
face. 

The  AEM  sounding  results  for  four  flights  over 


floe  site  1  are  shown  in  Figure  29.  The  average 
measured  S-I  thickness  of  lines  C  through  I  (C-I), 
a  swath  width  of  30  m,  is  also  shown.  It  is  apparent 
that  the  AEM  thickness  profiles  follow  the  meas¬ 
ured  long-period  S-I  thickness  trends  but  not  the 
short-period  undulations.  As  previously  discussed, 
this  is  due  to  the  averaging  associated  with  the 
footprint  area  over  which  the  ice  thickness  is  being 
integrated  as  well  as  the  use  of  a  one-dimensional 
model  to  characterize  three-dimensional  relief.  At 
present,  short-period  thickness  variations,  which 
occur  over  distances  of  less  than  two  or  three  times 
the  bird  height  above  the  surface,  are  not  well  de¬ 
fined,  if  at  all,  in  the  AEM  sounding  profiles. 

Sea  ice  thickness  requirements  vary  depending 
on  the  user's  need.  A  very  local  site-specific  ice 
thickness  is  preferred  by  someone  wanting  thin  ice 
for  instrument  deployment  or  a  thick  ice  site  to 
support  a  heavy  object.  An  average  of  the  relative 
sea  ice  thickness  variation  over  a  very  large  area 
would  be  extremely  useful  to  someone  interested 
in  modeling  ice  pack  movement  dynamics  or  in  as¬ 
sessing  heat  transfer  from  an  ice-covered  ocean  to 
the  atmosphere.  AEM  sounding  would  appear  to 
provide  S-I  thickness  information  well  suited  to 
the  latter  user  needs. 

All  AEM  profiles  appear  to  track  the  measured 
S-1  thickness  relief  along  the  entire  grid  length,  ex¬ 
cept  for  those  of  flight  FI  2L3.  This  profile  indicates 
greater  thicknesses  beyond  a  grid  distance  of  about 
310  m.  System  drift  was  again  the  problem. 

Similar  presentations  of  the  flight  profile  data 
taken  over  floe  sites  2  and  3  are  presented  in  Fig¬ 
ures  30  and  31 .  Again,  note  the  smoothing  of  the  ice 
relief  in  the  AEM  sounding  results.  Drift  also  ap¬ 
pears  to  have  affected  the  profile  results  over  floe 
3.  In  Figure  31 ,  flight  line  FI  3L5 's  ice  thickness  pro¬ 
file  is  seen  to  be  higher  than  the  average  measured 
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Figure  29.  Ai’crn^c  snow  plus  icc  thickness  for  30-in-widc  swath  C-I  down  center  of;^rid  on  MY  floe  I  ns  AEM-determined 
snow  plus  ice  thickness  for  four  sounding  fl>S^ds  (F...). 
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Figure  30.  AiK'ragc  sitoiv  plus  ice  thickness  for  30-m-unde  sivath  C-I  down  center  of  grid  on  MY  floe  2  vs  AEM- 
determined  snoiv  plus  ice  thickness  for  four  sounding  flights  (F...). 
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Figure  31.  Average  snoiv  plus  ice  thickness  for  30-nt-ivide  swath  C-I  down  center  of  grid  on  MY  floe  3  vs  AEM- 
determined  snoiv  plus  ice  thickness  for  three  sounding  flights  (F...). 
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Figure  32.  Average  snow  plus  ice  thickness  for  30-in-ivide  swath  C-I  down  center  of  grid  at  the  FY  ridge  vs  AEM-deterinined 
snow  plus  ice  thickness  for  three  sounding  flights  (F...). 
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values  (the  C-I  profile  line)  at  the  beginning  of  the 
grid,  while  those  for  flight  line  F13L6  and  F13L8 
are  higher  than  the  C-I  profile  line  at  the  end  of  the 
grid. 

The  AEM  profiles  taken  at  the  first-year  pressure 
ridge  site  are  shown  in  Figure  32.  Tlie  most  striking 
aspect  of  the  AEM  profiles  is  that  they  do  not  ade¬ 
quately  define  the  thick  ice  of  the  ridge.  Here  too, 
because  of  footprint  size  and  model  constraints,  as 
well  as  the  effect  of  the  seawater-filled  submerged 
ice  block  structure  and  the  relatively  steep-sided 
blocky  relief  of  the  ridge  keel,  it  is  currently  not 
possible  to  properly  define  such  short-period  relief 
that  occurs  over  a  distance  of  two  to  three  times  the 
bird  height  above  the  surface.  Since  the  thick  ridge 
ice  area  is  less  than  60  m  wide  and  the  bird  was 
flown  at  about  23  m  above  the  undeformed  surface 
relief,  the  results  shown  in  Figure  32  are  not  un¬ 
expected.  The  AEM  profiles  do  agree  with  the 
measured  thickness  of  the  thinner  ice  on  each  side 
of  the  ridge  and  do  indicate  the  location  of  the 
thicker  ridge  ice.  Smoothing  of  the  relief  in  the  area 
of  the  ridge  is  again  associated  with  footprint  size, 
ridge  geometry,  seawater  between  the  ice  blocks 
of  the  keel  and  1-D  model  constraints. 

Several  methods  were  used  to  evaluate  how 
well  the  AEM  S-I  thickness  determinations  agreed 
with  the  drill  hole  measurements.  Information 
related  to  each  AEM  sounding  flight  at  each  grid 


Table  10.  Average  AEM  S-I  thickness  and  related 
standard  error  for  each  sounding  flight  down  grids 
on  multiyear  floes  1,  2,  and  3  and  at  the  first-year 
ridge  site. 


Grid 

i^itc 

Grid 

length 

linilcf) 

Flight 

Flight  Speed 
no.  (km/hr) 

No.  of 
^^mnding^ 

Mean  AEM 
S-I  thickness 
(tn) 

Std. 

del’. 

On) 

Flex,'  1 

450 

F12L3 

77 

22 

4.27 

1.09 

F13L0 

52 

32 

4.29 

0.86 

F13L1 

54 

31 

4.09 

0.92 

F14L7 

43 

37 

4.11 

0.88 

Roe  2 

200 

F12L5 

60 

13 

5.84 

0.99 

F13L2 

51 

15 

5.64 

0.75 

F13L3 

60 

13 

5.68 

0.91 

F13L4 

42 

17 

5.76 

0,87 

Rexj  3 

200 

F13L5 

36 

21 

3.65 

1.33 

F13L6 

42 

18 

4.30 

1.23 

F13L8 

45 

17 

4.51 

1.14 

FYridg 

e  230 

FllLO 

67 

13 

3.61 

2.69 

F14L1 

25 

33 

3.25 

2.16 

F14L3 

21 

40 

2.74 

2.31 

site  is  listed  in  Table  1 0,  and  a  comparison  between 
the  average  drill-hole-measured  S-I  thickness  and 
the  average  AEM-determined  S-I  thickness  for 
each  sounding  flight  is  listed  in  Table  11. 

The  results  in  Tables  lOand  11  for  floes  1,2,  and 
3  do  not  provide  a  clear  indication  of  the  effect  of 
flight  speed,  and  therefore  the  number  of  sound¬ 
ings,  on  the  variation  between  the  average  AEM- 
determined  S-I  thickness  and  the  measured  value. 
One  would  expect  that  better  agreement  would 
have  occurred  as  more  AEM  soundings  were  ob¬ 
tained  along  a  grid.  Indeed,  just  the  opposite  is  ap¬ 
parent  for  the  flights  down  floe  3  and  at  the  FY 
ridge  site.  However,  these  may  have  been  the  re¬ 
sult  of  one  or  more  inconsistencies.  For  example, 
one  flight  line  did  not  exactly  follow  the  track  of 
another,  the  flight  elevations  and  therefore  foot¬ 
print  size  varied,  and  because  electronic  drift  was 
not  necessarily  linear  for  each  sounding  run,  the  S- 
I  thickness  determinations  may  be  slightly  differ¬ 
ent.  The  effect  of  flight  speed  on  the  sounding  re¬ 
sults  therefore  needs  further  study. 

The  data  in  Table  1 1  do  indicate  that  for  ice  floes 
with  low  to  moderate  relief,  AEM  sounding  is  cap¬ 
able  of  providing  a  good  assessment  of  an  ice  floe's 
average  S-I  thickness.  Indeed,  it  seems  that  once 
drift  and  calibration  problems  are  resolved  it  should 
be  possible  to  determine  the  average  S-I  thickness 
of  ice  floes,  with  "moderate"  relief,  to  within  5%. 

Table  11.  Percent  difference  between  average 
measured  S-I  thickness  for  lines  C-I  vs  average 
S-I  thickness  determined  by  AEM  sounding  for 
each  flight  down  the  grids  on  multiyear  floes  1, 2, 
and  3  and  at  the  first-year  ridge  sites. 


Grid 

site 

no. 

Measured 
S-I  thickness 
On) 

AEM 

S-I  thickness 
On) 

Difference 
Meas  i’s  AEM 
(’;)) 

Rex;  1 

F12L3 

3.93 

4.27 

+9 

F13L0 

4.29 

+9 

F13L1 

4.09 

+4 

F14L7 

4.11 

+5 

Rex;  2 

F12L5 

5.64 

5.84 

+4 

F13L2 

5.64 

= 

F13L3 

5.68 

- 

F13L4 

5.76 

+2 

Rex;  3 

F13L5 

4.22 

3.65 

-14 

F13L6 

4.30 

+2 

F13L8 

4.50 

+7 

FY  ridge' 

FllLO 

3.94 

3.61 

-8 

F14L1 

3.25 

-18 

F14L3 

2.74 

-30 
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Another  assessment  between  the  measured  and 
AEM  determined  S-I  thickness  was  made  by  com¬ 
paring  the  measured  vs  AEM  thickness  values  at 
stations  spaced  25  m  apart  along  the  center  grid 
line  on  floes  1  and  2.  Hoe  1  was  chosen  because  of 
its  relatively  low-lying  relief  and  floe  2  because  of 
its  thicker,  more  complex  relief.  The  S-I  thickness 
at  each  station  was  determined  by  averaging  the 
drill  hole  measurements  along  lines  C  through  I  (a 
30-m-wide  spacing)  and  for  the  row  on  which  the 
station  occurred  as  well  as  on  the  three  rows  before 
and  after  the  station  (a  30-m-wide  spacing).  Thus, 
the  S-I  thickness  used  for  each  station  was  the 
average  of  the  drill  hole  measurements  made  on  a 
5-m  grid  within  a  30-m-square  area  centered  at  the 
station.  This  average  or  measured  S-I  thickness  for 
each  station  was  used  because  it  is  a  reasonable  ap¬ 
proximation  to  the  integrated  thickness  sounded 
within  the  AEM  footprint.  Since  there  was  no  ice 
thickness  information  taken  before  the  first  or 
after  the  last  station  on  each  grid,  these  two  stations 
were  not  included  in  the  assessment. 

The  measured  and  AEM-determinedstation  S-I 
thickness  for  each  flight  over  the  grid  on  floes  1 
and  2  are  listed  in  Tables  12  and  13,  respectively. 


Also  listed  is  AT,  the  difference  of  the  AEM-de- 
termined  S-I  thickness  from  the  measured  value. 
The  most  significant  difference  occurred  at  the  25- 
m  station  on  flight  F13L0  over  floe  1.  Here  the 
AEM-determined  thickness  value  was  37%  above 
the  measured  one.  However,  at  most  stations  the 
difference  is  less  than  15%  .  The  standard  deviation 
of  the  AT  values  for  each  flight  is  also  listed.  These 
statistical  values  are  an  indication  of  the  relative 
dispersion  that  occurred  between  the  measured 
and  the  AEM-determined  i  -1  thicknesses.  Tire  dis¬ 
persion  is  numerically  larger  for  the  flights  over 
the  rougher  relief  of  floe  2,  but  in  relation  to  indi¬ 
vidual  station  S-I  thickness  the  deviations  are  quite 
similar,  again  about  15%'.  In  short,  the  AEM  S-I 
thickness  determination  for  each  grid  station  was 
generally  within  15%.  of  the  measured  value. 

Fhght  line  F15L1  ran  from  near  Gull  Island  to 
the  southeast  shore  of  Prudhoe  Bay  (Fig.  1  and  6). 
A  portion  of  this  line  included  ice  that  had  grown 
down  and  became  frozen  to  the  shallow  seabed. 
The  ice  cover  along  the  line  was  not  deformed  with 
pressure  ridges.  The  AEM-determined  S-I  thick¬ 
ness  and  sub-ice  water  depth  along  the  flight  line 
are  shown  in  Figure  33.  Large  variations  are  clearly 


Table  12.  Measured  and  AEM  determined  S-I  thickness  for  the  25-m  stations  along 
the  center  line  of  the  grid  on  multiyear  floe  1. 


AEM  flight 


Measured 

F12L3 

F13L0 

FnU 

FUL7 

Dist.  S-f  thickness 

S-I  thickness 

AT* 

S-I  thickness 

AT 

S-I  thickness 

AT 

S-I  thickness 

AT 

(m)  (m) 

(in) 

(in) 

(ill) 

(in) 

(m) 

(III) 

(in) 

(in) 

0 


25 

3.68 

4.80 

+1.12 

5.04 

+1.36 

4.39 

+0.71 

4.54 

+0.86 

50 

4.27 

3.94 

-0.33 

4.36 

+0.09 

4.13 

-0.14 

4.55 

+0.28 

75 

3.79 

3.11 

-0.68 

3.88 

+0.09 

4.13 

+0.34 

4.04 

+0.25 

100 

2.49 

3.06 

+0.57 

2.96 

+0.47 

3.02 

+0.53 

3.22 

+0.73 

125 

3.11 

3.06 

-0.05 

3.24 

+0.13 

2.86 

-0.25 

3.09 

-0.02 

150 

2.73 

2.91 

+0.18 

3.31 

+0.58 

3.02 

+0.29 

3.06 

+0.33 

175 

2.72 

2.79 

+0.07 

3.24 

+0.52 

3.02 

+0.30 

2.93 

+0.21 

200 

2.73 

2.68 

-0.04 

3.28 

+0.55 

2.89 

+0.16 

3.24 

+0.51 

225 

3.43 

3.98 

+0.55 

3.71 

+0.28 

3.89 

+0.46 

3.48 

+0.05 

250 

4.77 

4.75 

-0.02 

4.77 

0.00 

4.66 

-0.11 

4.32 

-0.45 

275 

4.98 

4.98 

0.00 

5.31 

+0.33 

4.66 

-0.32 

4.87 

-0.1 

300 

4.05 

4.25 

+0.20 

4.48 

+0.43 

4.00 

-0.05 

4.45 

+0.40 

325 

3.69 

4.25 

+0.56 

4.10 

+0.41 

4.00 

+0.31 

4.18 

+0.49 

350 

3.70 

4.89 

+  1.19 

475 

+0.55 

4.42 

+0.72 

4.18 

+0.48 

375 

4.73 

5.29 

+0.56 

4.72 

-0.01 

4.42 

-0.31 

4.61 

-0.12 

400 

4.71 

5.15 

+0.44 

5.34 

+0.63 

4.60 

+0.63 

5.07 

+0.36 

425 

450 

4.69 

5.85 

+1.16 

4.81 

+0.12 

5.09 

+0.40 

4.77 

+0.08 

Mean 

3.78 

4.10 

+0.32 

4.16 

+0.38 

3.94 

+0.17 

4.04 

+0.25 

Std.dev. 

0.83 

1.01 

0.52 

0.78 

0.33 

0.72 

0.34 

0.71 

0.33 

*  AT  is  difference  between  measured  and  AEM  determined  S-I  thickness  for  related  flight  (F . ). 
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Table  13.  Measured  and  AEM  determined  S-I  thickness  for  the  25-m  stations  along 
the  center  line  of  the  grid  on  multiyear  floe  2. 


_ AEM  fliaht _ 

Mciisiin-d  F12L5  F12L2  FT3L3  F13L4 


Dist. 

S-I  thickness 

S-I thickness 

AT* 

S-I  thickness 

AT 

S-I  thickness 

AT 

S-I  thickness 

AT 

(m) 

(in) 

(in) 

(in) 

<m) 

(in) 

(III) 

(in) 

(in) 

(in) 

0 


25 

4.75 

5.11 

+0.36 

4.62 

-0.13 

5.11 

+0,36 

5.05 

+0.30 

50 

6.39 

6.88 

+0.29 

5.57 

-0.82 

7.31 

0.92 

6.53 

+0.14 

75 

7.79 

6.88 

-0.91 

7.43 

0.36 

6.44 

-1.35 

6.67 

-1.12 

UX) 

-5.64 

6.08 

+0.44 

6.08 

+0,44 

5.80 

+0.16 

6.67 

+1.03 

125 

5.07 

5.49 

_0.42 

5.48 

.0.41 

5.53 

+0.46 

5.38 

+0.31 

150 

5.27 

5.49 

+0.22 

5.48 

0.21 

5.14 

-0.13 

6.39 

+1.12 

175 

200 

5.34 

5.01 

-0.33 

5.48 

+0.14 

5.30 

-0.04 

5.83 

+0.49 

Mean 

5.75 

5.85 

+0.07 

5.73 

-0.02 

5.80 

-0.21 

6.07 

+0.32 

Std.  dev. 

1.04 

0.78 

0.51 

0.86 

0.45 

0.81 

0.68 

0.66 

0.74 

•  AT  is  difference  between  measured  and  AEM  determined  S-I  thickness  for  relatc'd  flight  (F . ). 


seen  in  the  S-I  thickness  profile  which  cannot  be 
explained  by  snow  drift  features  or  seabed  topog¬ 
raphy.  The  former  were  short-period  surface  relief 
features  of  less  than  0.5-m-height  and  the  latter  is 
a  mild  relief  feature  sloping  into  deeper  water.  The 
variations  shown  are  an  example  of  AEM  system 
noise  and  instability  effects,  which  occurred  unex¬ 
pectedly  or  were  not  adequately  attended  to  prior 
to  the  flight.  While  the  average  AEM-determined 
S-I  thickness  may  be  representative  of  that  along 
the  flight  line,  it  is  clear  that  improvements  in  sys¬ 
tem  performance  are  in  order. 


Flight  line  F12L1  was  a  flight  over  the  fast  ice, 
which  began  south  of  Reindeer  Island  and  ended 
near  the  Prudhoe  Bay  West  Dock  (Fig.  4).  Tlie  AEM 
sounding  data  for  the  southern  half  of  this  flight 
line  were  processed  to  provide  both  S-I  thickness 
and  sub-ice  water  depth  as  shown  in  Figure  34.  In 
the  analysis  of  the  AEM  data  the  conductivity  of 
the  seawater  was  set  at  2.6  S/m.  As  seen  in  Figure 
34,  stations  D,  E  and  F  were  over-flown.  Also 
shown  is  that  the  AEM-determined  S-I  thickness  is 
in  good  agreement  with  the  measured  value  at  sta¬ 
tions  E  and  F  and  that  the  AEM-determined  sub- 
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Figure  33.  AEM-dcterviiucd  snoiv  plus  ice  thickness  and  sub-ice  water  depth  for  flight  line  F15L4.  See 
Figures  1  and  4  for  line  location. 
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Fi^’i/rc’  34.  AEM-deicrmincd  mow  plus  ice  thickness  and  sub-ice  water  depth  and  a  composite  cross  section  of  these  measurements  for 
flight  line  F22U.  See  Figure  2  for  line  location. 


ice  water  depths  at  each  station  are  within  10%  of 
the  measured  ones.  There  is  a  significant  difference 
in  theAEM  vs  measured  S-I  thickness  at  station  D. 
This  could  be  the  result  of  system  drift,  which  may 
have  affected  the  data  after  about  fiducial  number 
4025  and  could  bethe  reason  for  the  anomalousS- 
I  thickness  area  between  fiducial  numbers  4110 
and  4120. 

Flight  line  FllLlO  (Fig.  6)  was  flown  over  a 
first-year  sea  ice  rubble  area  and  then  down  an 
approximately  700-m-long  section  of  20  cm  ± 2  cm 
thick  snow-free  lead  ice.  The  AEM  S-I  thickness 
profile  along  this  flight  line  is  shown  in  Figure  35. 
The  lead  began  at  about  fiducial  number  4850.  The 
average  ATM-determined  lead  ice  thickness  was 
23  cm,  which  is  in  good  agreement  with  the  meas¬ 
ured  value  given  above. 

Flight  line  F15L4,  as  located  in  Figure  6,  was  an 
approximately  20-km-long  sounding  run  made 
over  the  near-shore  sea  ice.  This  run  was  intended 


to  represent  the  type  of  AEM  sounding  flight  that 
would  occur  during  a  pack  ice  sounding  mission. 
The  S-1  thickness  profile  along  this  flight  line  is 
shown  in  Figure  36.  At  abou  t  fiducial  numbers 
4670  and  5150  there  were  no  S-I  thickness  data 
available  du^'  to  saturation  of  the  receivers  as  a  re¬ 
sult  of  flying  the  bird  too  low  (close  to  the  sea  wa¬ 
ter).  The  flight  line  was  limited  in  length  by  an 
offshore  cloud  bank  that  prevented  VFR  flying  far¬ 
ther  north. 

The  relative  distribution  of  the  S-I  ice  thickness¬ 
es  sounded  along  the  flight  line  is  shown  in  Figure 
37.  A  thin  lead  ice  component  is  revealed  in  the 
data,  as  is  the  smaller  amount  of  thick  pressure 
ridge  ice.  The  average  S-I  thickness  was  3.7  m. 
There  are  no  ground  truth  data  to  validate  the  S- 1 
thickness  distribution  shown.  In  addition,  the  flight 
line  probably  did  not  cross  ice  representative  of 
that  located  in  other  Arctic  Ocean  areas.  However, 
it  may  be  of  interest  to  note  that  Wittmann  and 
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Fif’iirc  37.  Histo^raiu,  probability  and  proportion  dis¬ 
tribution  of  snoio  plus  icc  thickness  along  flight  line 
F15L4. 


Schule  (1966)  estimated,  from  sonar  profiles  of 
under-ice  relief  obtained  during  submarine  cruises, 
that  the  average  arctic  pack  ice  thickness  is  about 
3.9  m,  which  is  in  agreement  with  the  more  recent 
finding  of  Colony  (1988)  of  3.8  m  obtained  from 
submarine  ice  bottom  profile  data  collected  in  the 
Beaufort  Sea.  Both  of  these  averages  are  in  agree¬ 
ment  with  the  average  3.7-m  AEM  value  given 
above.  In  any  event  the  trend  of  the  distribution  in 
Figure  37  is  as  expected  and  is  in  good  agreement 
with  the  many  reports  that  give  the  distribution  of 
ice  drafts  obtained  during  submarine  transects  un¬ 
der  the  Arctic  pack  ice  in  the  (e.g.,  Blidberg  et  al. 
1979  and  McLaren  1988). 

CONCLUDING  REMARKS 

'T'he  most  serious  EM  system  problem  encoun¬ 
tered  during  the  1987  field  work  was  drift,  particu¬ 
larly  in  the  high-frequency  channels.  This  drift 
was  mainly  caused  by  temperature-related  changes 
in  the  analog  receiver  console,  particularly  in  the 
band-pass  filter  assemblies.  Our  solution  to  this 
problem  will  be  to  use  a  digital  receiver  specifically 
designed  to  minimize  temperature  drift-related 
problems. 

A  secondary  source  of  drift,  and  one  which  is 
harder  to  pin  down,  was  temperature-related  drift 
in  the  bird's  transmitter  and  receiver  assemblies. 
As  transmitter  coils  warm  up  and  expand,  there  is 
a  slight  shift  in  transmitter  coil  frequency  and  mo¬ 
ment.  Frequency  drift  is  the  more  important  of 
these  two  changes.  The  receiver  coils  have  a  rela¬ 
tively  sharp  resonance,  with  a  Q  factor  of  10. 
Therefore,  a  small  frequency'  shift  in  transmitter 
frequency  will  result  m  a  significant  amplitude 
and  phase  change  in  the  received  signal,  which  in 
turn  leads  to  systematic  errors  in  the  ice  thickiress 
and  conductivity  estimates  as  well  as  in  the  bathy¬ 
metric  results.  Moment  drift  results  in  smaller  ab¬ 
solute  secondary  field  pickup  at  the  receiver,  which 
in  turn  causes  small  changes  in  both  phase  compo¬ 
nents.  Both  of  these  problems  must  be  resolved  in 
the  development  of  a  new  transmitter  module. 

An  annoying  but  less  important  problem  was 
the  presence  of  a  significant  amount  of  noise  in  the 
high  frequency  channel.  Tliis  noise  was  correlated 
with  motion  of  the  tow  cable  with  respect  to  the 
bird,  and  was  reduced  by  immobilizing  the  cables 
as  much  as  possible  near  the  bird  and  by  tying 
spoiler  flags  to  the  tow  cable  to  reduce  the  ampli¬ 
tude  of  wind-driven  oscillation.  However,  the  key 
source  of  noise  was  discovv  red  after  the  field  sur¬ 
vey.  A  few  inches  of  unshielded  wiring  was  found 
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in  the  bird.  Tliis  wiring  acted  as  an  electrostatic  an¬ 
tenna  that  picked  up  EM  radiation  noise.  Such 
shielding  errors  should  have  been  detected  and 
corrected  by  the  manufacturer,  had  proper  testing 
occurred  prior  to  shipping  the  AEM  system  to 
Alaska. 

A  most  important  problem  involved  the  behav¬ 
ior  of  the  laser  altimeter's  control  console.  As  men¬ 
tioned  earlier,  the  console  or  the  interface  between 
it  and  the  data  recording  system  delayed  logging 
of  laser  altitudedata  by2.5seconds.  Acombination 
of  effects  in  the  laser  console  and  the  data  acquisi¬ 
tion  system's  control  program  is  probably  respon¬ 
sible  and  will  be  resolved  in  the  future. 

Tire  problems  encountered  in  field-  and  post¬ 
processing  of  the  AEM  data  derived  almost  entirely 
from  a  single  source:  the  altimeter  time  delay.  Dur¬ 
ing  the  field  processing  effort,  the  systematic  incon¬ 
sistencies  in  the  data  set  caused  by  the  altimeter 
delay  led  to  many  fruitless  hours  of  searching 
through  the  entire  software  system  and  the  AEM 
data  lookingforpathologicalnonlineari  ties  in  drift. 
These  efforts  continued  until  the  altimeter  delay 
was  discovered.  Once  this  problem  was  dealt  with, 
it  was  straightforward  to  obtain  an  acceptable  sys¬ 
tem  calibration  based  on  ground  truth  (a  known 
snow-ice  thickness  and  water  depth). 

Tire  radar  altimeter  used  to  determine  helicopter 
elevation  would  frequently  stopdisplaying  altitude 
during  a  flight.  The  result  could  have  been  a  refusal 
by  the  pilot  to  fly  at  the  optimal  (bird)  survey 
height  of  20  m.  However,  he  was  willing  to  fly  us¬ 
ing  the  bird  altitudeasdetermined  by  the  laser  alti¬ 
meter.  It  seems  most  likely  that  the  radar  problem 
was  produced  by  a  combination  of  effects,  particu¬ 
larly  weak  or  diffuse  reflections  from  the  snow 
surface  and  the  lack  of  a  good  ground  plane  around 
the  backing  plate  of  the  radar  antenna  (the  ASTAR 
helicopter  used  had  a  composite  shell  with  a  low 
electrical  conductivity). 

The  failure  of  the  RMS  data  recording  system 
has  already  been  mentioned,  as  were  the  steps 
used  to  correct  it. 

The  main  thrust  of  this  airborne  sea  ice  meas¬ 
urement  program  was  to  design,  construct  and 
test  a  prototype  bird  having  a  3-m  nominal  coil 
separation  and  frequencies  ranging  from  8(X)  Hz  to 
50  kHz.  The  prototype  AEM  bird  was  very  stable 
in  flight  but  the  electronics  did  not  perform  as  ex¬ 
pected  during  arctic  testing.  With  the  corrections 
made  after  the  field  trials  and  modifications  pro¬ 
posed,  the  AEM  system  should  work  very  satisfac¬ 
torily. 

One  aspect  of  this  work  was  to  implement  and 


test  a  practical,  robust  and  fast  field-processing 
system  for  analyst  s  of  the  airborne  ice  measurement 
data.  Turnaround  times  of  less  than  one  day  were 
demonstrated  in  the  field,  and  have  been  further 
shortened  by  subsequent  software  innovations. 
The  quality  of  the  in-field  interpretations  did  not 
meet  expectations,  but  this  was  the  result  of  an  in¬ 
terfacing  problem  with  the  laser  altimeter,  not  of 
inadequacy  in  the  processing  software.  Certainly, 
the  final  interpretation  conforms  well  with  the 
ground-truth  measurements  wherever  system  drift 
was  not  an  undue  problem. 

Real-time  processing  was  one  of  the  ultimate 
objectives  of  this  program.  Designing  software  to 
perform  this  task  is  now  a  fairly  straightforward 
task,  which  mainly  requires  a  substantial  input  of 
time  and  effort.  However,  the  system  will  need  to 
have  sufficient  "intelligence"  to  cope  with  changes 
in  the  character  of  the  surface  below,  noise  due  to 
radio  transmissions  from  the  aircraft  and  occasional 
saturations  of  the  AEM  system  if  the  bird  is  flown 
too  low,  all  in  real  time.  If  drift  cannot  be  effectively 
eliminated  from  the  AEM  system,  then  thesoftware 
will  need  to  predict  and  remove  drift  effects  based 
on  previous  drift  behavior.  This  may  well  prove 
rather  difficult. 

A  number  of  other  objectives  were  also  achieved . 
For  example,  this  work  produced  estimates  of  sea 
ice  conductivity  using  an  airborne  measurement 
system.  While  these  results  were  erratic,  due  to 
system  drift  and  noise  problems,  they  certainly 
constitute  an  exciting  step  toward  airborne  sea  ice 
characterization. 

More  detailed  measurement  of  ice  conductivity 
over  a  wide  range  of  ice  types  will  require  some 
additional  instrument  development  effort,  because 
frequencies  of  200  kHz-300  kHz  or  higher  may  be 
necessary  (Becker  et  al.  1987).  A  100-kHz  subsystem 
has  now  been  installed  in  the  small  bird  but  has  not 
been  test-flown. 

One  approach  toward  development  of  an  inter¬ 
nal  precision  calibration  system  was  investigated 
and  found  to  be  inadequate.  Further  development 
in  this  area  is  considered  a  requirement  that  must 
be  satisfactorily  implemented  in  an  AEM  system 
before  the  teclmology  can  be  considered  for  routine 
sea  ice  tliickness  sounding  surveys.  Such  equip¬ 
ment  is  essential  for  a  production  system,  since 
ground-truth  calibrations  after  bird  modification 
or  maintenance  may  not  always  be  possible. 

It  is  clear  that  one-dimensional  interpretation 
of  sea  ice  thickness  works  well  for  relatively  flat 
ice,  but  when  major  changes  in  ice  thickness  occur 
over  lateral  distances  of  less  than  about  two  bird 
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elevations,  errors  will  be  introduced  by  the  1-D  ap¬ 
proximation  due  to  footprint  area  effects.  Thus, 
the  measurement  of  sea  ice  thickness  using  AEM 
techniques  may  no  longer  be  limited  by  instru¬ 
mentation  problems,  but  rather  by  inadequate  in- 
terpretational  capability  in  the  vicinity  of  compli¬ 
cated  2-D  or  3-D  ice  structures. 

Becker  et  al.  (1987)  and  Becker  and  Liu  (1988) 
have  studied  the  ice  measurement  problem  and 
identified  a  method  that  might  be  suitable  for  2-D 
or  3-D  ice  thickness  determination,  although  it 
does  not  appear  to  address  the  interpretation  of  ice 
conductivity  in  such  situations.  Implementing  this 
method  in  a  real-time  processing  system  may  be 
difficult,  since  it  requires  analysis  of  multiple  mea¬ 
surement  points  simultaneously  and  is  currently 
too  processing-time-intensive,  requiring  tens  of 
seconds  per  data  point  on  a  mainframe  computer. 

Even  if  a  fully  developed,  real-time  AEM  pro¬ 
cessing  algorithm  for  sea  ice  thickness  were  avail¬ 
able  for  estimating  the  cross-sectional  variation  of 
three-dimensional  sea  ice  structures,  it  would  not 
be  possible  to  fully  characterize  all  ice  morphology. 
This  would  be  especially  true  for  first-year  pressure 
ridges  and  rubble  formations  with  a  width  of 
about  one  bird  elevation  and  for  ice  formations 
with  complex  geometry.  As  shown  in  Figure  38, 
first-year  pressure  ridges  can  have  a  multitude  of 
cross-sectional  and  internal  geometries.  These 
structure  variations  give  rise  to  complex  secondary 
electromagnetic  field  responses  that  would  be  dif- 
ficult  to  interpret.  While  this  structure  variation 
certainly  becomes  less  severe  after  one  or  more 
melt  seasons,  it  is  not  realistic  to  expect  that  the 
keel  relief  depicted  in  Figure  38b,  38c  or  39  could 
be  well-defined.  Smoothing  of  the  compi?'''  ice  re¬ 
lief  will  occur,  with  the  deepest  structure  in  an  ice 
formation  being  poorly  determined.  Nevertheless, 
the  results  of  this  program  indicate  that  it  should 
be  possible  to  determine  the  mean  thickness  of 
arctic  pack  ice  to  wi  thin  10%.  Such  AEM  sounding 
capability  should  be  of  considerable  interest  to 
those  studying  the  thermal  and  motion  dynamics 
of  the  arctic  pack  ice  and  to  those  interested  in  lo¬ 
cating  thick  or  thin  sea  ice. 

This  study  and  our  previous  one  (Kovacs  et  al. 
1987a)  revealed  that  an  AEM  sounding  system  is 
capable  of  measuring  undeformed  sea  ice  thickness 
from  a  few  centimeters  to  5  or  more  meters  thick  to 
within  10%  of  the  drill  hole  measured  value.  Accu¬ 
racies  within  5%  can  be  expected  when  system  ge¬ 
nerated  noise  is  eliminated  and  drift  is  properly 
controlled.  However,  conventional  AEM  systems 
use  analog  receivers  which  tend  to  drift  with  time 


Figure  38.  Examples  of  idealized  first-year 
pressure  ridge  cross  sections  and  their  trans¬ 
formation  into  multiyear  features.  Cross  sec¬ 
tions  a,  band  care  modified  from  Grishchenko 
(2988). 


and  temperature  variations.  The  result  is  that  pre¬ 
cise  estimates  of  the  secondary  magnetic  field  re¬ 
sponse  from  the  seawater  become  less  reliable  as 
the  time  of  the  ice  sounding  survey  increases.  In 
addition,  the  number  of  frequencies  (transnaitter- 
receiver  coil  pairs)  that  can  be  installed  in  a  small 
bird  is  limited.  This  limitation  prevents  selection 
of  another  frequency  during  a  survey  flight  should 
one  of  the  built-in  frequencies  become  affected  by 
some  interference  noise. 

To  overcome  these  problems,  we  are  developing 
a  system  with  a  single  transmitter  coil  that  will  ra¬ 
diate  a  broad  spectrum  of  electromagnetic  energy 
and  a  receiver  coil  that  will  be  used  to  measure  the 
secondary  magnetic  field  over  a  wide  frequency 
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Figure  39.  Inclined  sea  ice  structure  under  a  multiyear  pressure  ridge. 


range.  The  use  of  a  single  transmitter  should  elim¬ 
inate  the  cross-frequency  noise  that  occurs  in  con¬ 
ventional  multicoil  pair  AEM  systems.  In  addition, 
digital  signal  processing  will  be  used  in  the  new 
system  to  minimize  and  possibly  eliminateconsole 
drift.  The  console  will  have  built-in  computer  hard¬ 
ware,  capable  of  extremely  rapid  calculations, 
which  will  allow  inflight  data  processing  and  dis¬ 
play  of  S-I  thickness. 
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APPENDIX  A:  BLOWUP  OF  FIGURE  28 
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S-l,  Snow-ice  Thickness  (m) 
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